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ABSTRACT 

We present new identifications of infrared counterparts to the population 
of hard X-ray sources near the Galactic center detected by the Chandra X-ray 
Observatory. We have spectroscopically confirmed 16 new massive stellar coun- 
terparts to the X-ray population, including nitrogen-type (WN) and carbon-type 
(WC) Wolf-Rayet stars, and O supergiants. These discoveries increase the to- 
tal sample of massive stellar X-ray sources in the Galactic center region to 30 
(possibly 31). For the majority of these sources, the X-ray photometry is con- 
sistent with thermal emission from plasma having temperatures in the range of 
kT = 1-8 keV or non-thermal emission having power-law indices in the range 
of — 1 < r < 3, and X-ray luminosities in the range of Lx ~ 10^^- 10'^^ erg s~^ 
(0.5-8.0 keV). Several sources have exhibited X-ray variability of several factors 
between observations. These X-ray properties are not a ubiquitous feature of 
single massive stars but are typical of massive binaries, in which the high-energy 
emission is generated by the collision of supersonic winds, or by accretion onto a 
compact companion. However, without direct evidence for companions, the pos- 
sibility of intrinsic hard X-ray generation from single stars cannot be completely 
ruled out. The spectral energy distributions of these sources exhibit significant 
infrared excess, attributable to free-free emission from ionized stellar winds, sup- 
plemented by hot dust emission in the case of the WC stars. With the exception 
of one object located near the outer regions of the Quintuplet cluster, most of the 
new stars appear isolated or in loose associations. Seven hydrogen-rich WN and 
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O stars are concentrated near the Sagittarius B H ll region, while other similar 
stars and more highly evolved hydrogen-poor WN and WC stars lie scattered 
within ^ 50 pc, in projection, of Sagitarrius A West. We discuss various mecha- 
nisms capable of generating the observed X-rays and the implications these stars 
have for massive star formation in the Galaxy's Central Molecular Zone. 



1. Introduction 

The Galactic center harbors the largest concentration of X-ray sources in the sky. In- 
deed, a total of 9017 X-ray sources have been detected in deep observations of the central 
2° X 0?8 of the Galaxy with the Chandra X-ray Observatory (Muno et al. 2009). The major- 
ity of these sources are highly absorbed by interstellar gas and dust, indicative of distances 
near the Galactic center {D ^ 8 kpc; Reid 1993). A wide variety of stellar phenomena 
comprise the Galactic center X-ray population, including accreting magnetic white dwarfs 
in cataclysmic variables (CVs; e. g., see Muno et al. 2003) and symbiotic binaries; accreting 
neutron stars and black holes in low- and high-mass X-ray binaries (e. g., see Pfahl, Rappa- 
port & Podsiadlowski 2002); Wolf-Rayet and 0-supergiant binaries with colliding supersonic 
winds; and, perhaps, rogue compact objects accreting from the dense Galactic center medium 
(Agol & Kamionkowski 2002). All such objects arc rare among stellar populations, and in- 
volve exotic forms of matter and high-energy radiation mechanisms, of which our knowledge 
is fragmentary. Thus, the characterization of the Galactic center X-ray population is an im- 
portant step toward improving our understanding of the end stages of stellar evolution, and 
the 'zoology" of X-ray sources contained within the only galactic nucleus we are currently 
able to resolve in such detail. 

X-ray properties alone, however, are insufficient to characterize the Galactic center X- 
ray population. The detection of counterparts at other wavelengths is necessary, although 
such a pursuit is limited by the high extinction suffered by starlight traversing 8 kpc from the 
Galactic center to Earth (Ay ~ 30 mag). Fortunately, the brightness distribution of infrared 
counterparts can be used to constrain the relative contributions of low- and high-mass objects 
to the sample. For instance, CVs with cool dwarf donors in the Galactic center should have 
very faint IR counterparts with K > 22 mag, while main-sequence stars earlier than BOV, 
and Wolf-Rayet stars, will have 7^ < 15 and < 12 mag, respectively. It was demonstrated 
in Mauerhan et al. (2009; hereafter referred to as Paper 1) that only ^ 6%±2% of the 6760 
absorbed X-ray sources have real infrared counterparts with Kg < 15.6 mag. Although this 
result is consistent with earlier studies that concluded the population to be widely dominated 
by CVs (Muno et al. 2003, Laycock et al. 2005), it suggests, nonetheless, that 100-300 
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of the X-ray sources in our sample have real infrared counterparts within a brightness range 
consistent with that of the population of late-type giants and hot massive stars observable 
in the Galactic center. Indeed, several X-ray sources with hot supergiant counterparts have 
already been discovered, including O supergiants, and nitrogen-type (WN) and carbon-type 
(WC) Wolf-Rayet stars (Muno et al. 2006, Mikles et al. 2006, Mauerhan, Muno & Morris 
2007, Hyodo et al. 2008). 

The successful detection of massive-star counterparts to Chandra sources is not sur- 
prising, given the environmental conditions at the Galactic center. The region boasts the 
highest star formation rate density in the Milky Way (10"'^ Mq yr^^ pc~^; Figer et al. 2004), 
owing to the immense reservoir of molecular gas that occupies the Central Molecular Zone 
(e.g., see Morris & Serabyn 1996). This material is vigorously forming massive stars near 
the Sagitarrius B H ll region (e. g., see dePree et al. 1998), and within the last several Myr 
has produced at least three extraordinary stellar clusters that are among the most massive 
and dense in the Galaxy: the Arches and Quintuplet (Nagata et al. 1995; Cotera et al. 1996, 
Figer et al. 1999, 2002), and the Central parsec cluster (Krabbe et al. 1995). These clusters 
are rich in Wolf-Rayet stars and O supergiants. Such stars are typically sources of soft, 
thermal X-rays with kT < 1 keV, while only a small subset of massive stars (those in close 
binaries) will produce hard X-rays with kT > 1-2 keV. Since the softer X-ray photons are 
heavily absorbed by interstellar gas and dust, only the harder sources are hkely to be detected 
at the Galactic center. This explains why the Arches and Quintuplet clusters contain only 
several detectable X-ray sources among their hundred or so massive stars (Yusef-Zadeh et al. 
2002; Law & Yusef-Zadeh 2004; Wang et al. 2006). Outside of these clusters, the relatively 
isolated X-ray-emitting supergiants discovered by Mauerhan, Muno & Morris (2007), which 
reside within ^ 10 projected pc of the Arches and Quintuplet clusters, may have been dy- 
namically ejected from one of these systems. Alternatively, these and other isolated massive 
stars throughout the Galactic center may be the products of an isolated mode of massive star 
formation operating in the region, in tandem with the formation of stellar clusters. Thus, the 
identification of infrared counterparts to the Chandra X-ray population can provide insight 
into the evolution of stellar clusters, and highlight previously unknown regions of massive 
star formation in the Galactic center. 

In this work, we present the results of our pursuit of counterparts via infrared spec- 
troscopy. In Section 2, we list the various facilities and instrumentation used for this work 
and describe our observations. In Section 3, the spectra of the newly found massive stars 
are presented, including new spectra of several previously discovered objects. In Section 4, 
the stellar parameters of confirmed sources are examined via their photometry and infrared 
spectral energy distributions (SEDs). In Section 5, the X-ray photometric properties of the 
massive stars are listed and compared with the greater field population of X-ray sources in 
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the Galactic center. In Section 6, we discuss potential hard-X-ray emission mechanisms, 
consider the spatial distribution, origin, and formation mode of these stars, and discuss the 
implications for the greater population of massive stars in the Central Molecular Zone. 

2. Observations 

We performed infrared spectroscopy of candidate counterparts to X-ray sources that 
were selected via the cross-correlation of the Chandra catalog of Muno et al. (2009) and the 
SIRIUS and Two Micron All Sky Survey (2MASS; Cutri et al. 2003) near-infrared catalogs. 
The details of the sample and the cross-correlation analysis are presented in Paper I. We 
targeted sources that (1) have highly absorbed soft X-ray fluxes and reddened near-infrared 
photometry (the so-called red infrared matches to hard X-ray sources from Paper 1), which 
is the characteristic of objects lying at the distance to the Galactic center, (2) lie at angular 
distances >7' from Sgr A*, thus avoiding the high stellar confusion and large number of 
accidental infrared/X-ray matches that lie within this radius (3) have Kg < 12 mag, and 
(4) have excellent X-ray astrometry, with positional uncertainties <1". We have currently 
obtained spectra of 52 infrared matches to X-ray sources. We observed in the K-hand since 
it contains many diagnostic features of both hot and cool stars (e.g., see Hanson et al. 1996, 
2005; Morris et al. 1996; Figer, McLean & Najarro 1997), and suffers the least amount of 
extinction in the near infrared. Eighteen sources, including those from Mauerhan, Muno & 
Morris (2007), exhibit characteristics of massive stars (see Section 3), while the remainder 
exhibit CO band-head absorption features near A2.3 /xm, and sodium absorption features 
near A2.21 /xm, typical of cool, late- type stars (spectra not presented here). In general, we 
assumed that a late-type match to an X-ray source was a spatial coincidence, since late- 
type stars vastly outnumber hot stars in the field. Thus, late-type counterparts will not be 
presented in this paper. However, we note that symbiotic X-ray binaries consisting of late- 
type giants or supergiants with accreting compact companions could be indistinguishable 
from random coincidences if emission-line accretion signatures were too weak to be detected 
in our spectra. Thus, we cannot completely rule out the presence of symbiotic stars in the 
late-type sample. 

We do not claim that the true statistics of infrared counterparts to X-ray sources are 
indicated by the number of massive-star detections presented below; there are several reasons 
for this: our spectroscopy campaign began before the completion of the Chandra survey 
and the final source catalog, which in the end allowed us to refine the X-ray point-source 
positions and uncertainties, and eliminate spurious detections (the method of refining the 
X-ray astrometry is discussed in Muno et al. 2009). Thus, several matches for which we had 
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obtained prior spectroscopy did not remain in the final candidate counterpart list that was 
based on the refined X-ray catalog. In fact, only 15 of the 35 sources which we confirmed 
spectroscopically to be late-type stars remained in the final list. Furthermore, candidate 
counterparts were also given priority for spectroscopy based upon their location in particular 
regions of interest. For instance, if a candidate's position on the sky places it near an 
interesting diffuse structure, such as a bubble-like or shell-like feature in the Spitzer /YRKC 
A8 //m images of the region, then it was given priority over other candidate counterparts in 
the hst. The motivation for this additional selection criterion hes in the fact that intense 
winds from massive stellar X-ray sources can sweep up and heat interstellar material. Thus, 
we interpreted the spatial association of candidate infrared/X-ray counterparts with diffuse 
mid-infrared structures as indicative of a higher likelihood that the X-ray source is associated 
with a massive star. Although this technique probably increased our success rate, it almost 
certainly imposed a selection effect that hampers our ability to statistically analyze our 
spectroscopic success rate for finding counterparts to the Chandra population. 

Table 1 lists the X-ray sources and associated massive-star counterparts that we targeted 
for spectroscopy, along with the instrument used, and the dates of the observations. The 
table also includes information on seeing and sky conditions. All sources are presented 
with their full Chandra designation in Table 1, and are hereafter referred to in the text 
by an abbreviated version (e.g., CXOGC J174532. 7— 285617 will hereafter be referred to as 
X174532.7). 

2.1. NASA Infrared Telescope Facility (IRTF) and SpeX 

Spectra of four stars were obtained using the SpeX medium-resolution spectrograph 
on the 3 m Infrared Telescope Facility (IRTF) telescope (Rayner et al. 2003), located on 
the summit of Mauna Kea in Hawaii. For X174656.3 and X174711.4, SpeX was used in 
the high-throughput, low-resolution prism mode, which provides coverage of the entire JHK 
bandpass, although we will only include the X-band in this work. A slit width of 0'.'3 was 
used, providing a spectral resolution ol 250. For the infrared counterparts to X174645.2 
and X174550.6, medium-resolution spectra were obtained using the short cross-dispersed 
mode (SXD) of SpeX, with slit widths of O'.'S and 0'.'3 for the respective stars, providing 
spectral resolutions of i? ~ 1200 and R ^ 2000 in the K-hand. All spectra were acquired 
in an ABBA nodding sequence in order to subtract the sky background, and to suppress 
the contribution of bad pixels. The spectra were reduced and extracted using the IDL-based 
software package Spextool, specially designed for the reduction of data obtained with SpeX 
on the IRTF (Gushing et al. 2004). For telluric correction, spectra of the AOV standard 
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star HD 162220 were obtained and applied to all of the low-resolution data, while the AOV 

star HDl 55379 was used for the medium- resolution spectra. Telluric correction was executed 
using the IDL package xtellcor (Vacca et al. 2003), which applies and removes model H I 
absorption lines from the AOV standard star before application to the science data. 

2.2. The United Kingdom Infrared Telescope (UKIRT) 

The 3 m United Kingdom Infrared Telescope (UKIRT) on Mauna Kea was used to ob- 
tain spectra of eight infrared counterparts: X174519.1, X174532.7, X174537.3, X174628.2, 
X174703.1, X174712.2, X174713.0, and X174725.3. The UKIRT 1-5 /xm Imager Spectrom- 
eter (UlST; Ramsay et al. 2004) was used in service mode as part of the UKIRT Service 
Programme. The short-K grism and 4-pixel slit (0'.'48) were used, which provided a spec- 
tral resolution of i? ~ 2133 and a wavelength range of 2.01-2.26 /xm. Spectra were ac- 
quired in an ABBA nodding sequence for sky subtraction. Basic reductions were executed 
by the Starlink ORACDR pipeline and the spectra were extracted using the IRAF routine 
APALL. The telluric spectrum of HD162220 was applied to all UKIRT/UIST data using the 
xtellcor_general program, a generalized version of the same program used for IRTF SpeX 
data (used for all stars not observed at the IRTF). 

2.3. The Anglo-Australian Telescope and IRIS2 

Spectra of X174516.1, X174616.6, and X174617.7 were obtained using the 4.1 m Anglo- 
Austrahan Telescope (AAT) on Siding Spring Mountain (Mount Woorat) in New South 
Wales, Australia. The IRIS2 instrument (Tinney et al. 2004) provided a spectral resolution 
of R ^ 2400 in the 7^-band, using the 1" slit. The data images were reduced using the 
Starlink ORACDR pipeline. The spectra were extracted using the IRAF routine APALL. The 
telluric spectrum was obtained from observations of the FOV star BS6441. 

2.4. The Southern Observatory for Astrophysical Reseeirch cind OSIRIS 

Spectra of X174502.8, X174508.9 and X174516.7 were obtained at the Southern Ob- 
servatory for Astrophysical Research (SOAR), located on Cerro Pachon in Chile. The Ohio 
State Infrared Imager Spectrometer (OSIRIS; Depoy et al. 1993) was used in high-resolution 
longslit mode, which provided i?~3000 in the A'-band. Stellar spectra were acquired in a slit- 
scan sequence of five positions separated by 5" each. The five spectral images were median 
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combined with a suitable bad-pixel rejection algorithm to produce a master sky spectrum, 
which was subtracted from each individual exposure. The AOV standard HD 155379 was 
used for telluric calibration. 



3. Spectral Classification 

In the following subsections, we estimate the spectral types of the infrared counterparts 
to Galactic center X-ray sources. All of the spectral criteria we use to classify these stars were 
adopted from previous spectroscopic studies of massive stars in the near-infrared, conducted 
by Morris et al. (1996), Hanson et al. (1996, 2005), Figer, McLean & Najarro (1997), and 
Martins et al. (2008). The central wavelengths of all spectral lines referred to in the text 
are adopted from Figer, McLean & Najarro (1997, their Table 2). 

3.1. Late O Supergiants 

Figure 1 shows the K-hand spectra of the infrared counterparts to X174537.3, X174502.8, 
and the Quintuplet member X174616.6 (star qF344 of Figer, McLean & Morris 1999). Each 
of these stars exhibits absorption lines of Br7 at A2.1661 fim and He I at A2.058, A2.113, 
and A2.1647 fim, which arc typical features of late O and early B supergiants. In addition, 
each star also exhibits a weak emission feature of N ill, which appears at A2.115 /xm, on the 
red side of the He I absorption hue. The main feature distinguishing the stars of Figure 1 
is the varying depth of the Br7 line, which tends to increase with later spectral type. The 
Br7 feature is weakest for X174502.8, deepest for X174537.3, and is of intermediate strength 
for X174616.6. The latter star exhibits a weak absorption line of He ll as well, which in- 
dicates that it is somewhat hotter and of slightly earlier spectral type than X174537.3 and 
X174502.8. Based on a comparison of these spectra with the OB stars of Hanson et al. (1996, 
2005) , we classify X174502.8 and X174537.3 as supergiants within the O9I~B0I range, and 
X174616.6 as an 08-091 star. 

Wc note that the X-ray detection of X174502.8 is questionable. The source detection 
was flagged as spurious and not included in the catalog of Muno et al. (2009), but owing 
to the source's position near an interesting shell-like feature in Spitzer A8 /xm images, it 
was targeted for spectroscopy nonetheless. This implies that either some legitimate X-ray 
detections near the threshold of significance have been discarded by our selection algorithm, 
or the spatial correlation of the OB star with a truly spurious X-ray detection was a very 
unusual coincidence. 
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3.2. Ecirly O Supergiants 

The infrared counterparts of X174532.7, X174628.2, X174703.1, and X174725.3 all ex- 
hibit very similar i^-band spectra, presented in Figure 2. The spectra are dominated by a 
complex of blended emission lines near A2. 112-2. 115 /xm that is dominated by He I and N 
III, and might also may contain contributions from C ill and O ill. Blueward of this hne 
complex is weaker emission lines of C IV at A2.069 and A2.078 //m. Br7 is seen in emission 
in X174628.2, although it may be nebular feature, while there is marginal detection of Br7 
absorption in X174532.7. Both X174532.7 and X174703.1 exhibit He ll absorption at A2.189 
fim, which is either absent or below the noise level in the other two stars. All of these 
features are consistent with the i^'-band spectra of early O supergiants, specifically 04-61 
stars (e. g., see Martins et al. 2008) and we classify them as such. 

The Arches cluster of the Galactic center is rich with stars of 04-61 spectral type (Figer 
et al. 2002, Najarro et al. 2004, Martins et al. 2008). For comparison, accompanying the 
spectra in Figure 2 is a model spectrum for the 04-61 star F18 of the Arches cluster from 
Martins et al. (2008). This model was computed using the stellar atmosphere code CMFGEN 
(Hillier & Miller 1998). The model star has the following values for effective temperature, 
luminosity, mass-loss rate, and terminal wind velocity of these stars: Te// = 36900 K, 
L = IO^-^Lq, M = 3.2 X 10~^Mo yr~\ and v^o = 2150 km s~^ It is reasonable to assume 
that the stars in Figure 2 have similar properties. 

3.3. WNh and OIf+ Stars 

The A'-band spectra of X174656.3, Xr74711.4, and X174617.0 are presented in Figure 
3. The spectra of these stars are dominated by emission fines of Br7 and the A2. 112-2. 115 
//m complex of He I, N ill, C ill and O ill. The spectrum of X174617.0 was first presented in 
Mauerhan, Muno & Morris (2007; hereafter referred to as MMM07), where it was classified 
as an 06If"'" star, but we present a refined version of it here. The Br7 line of this particular 
star exhibits an asymmetry on its blue side owing to a contribution from He I, although 
the detectability of this feature may be due to the fact that the spectrum of X174617.0 
has significantly higher spectral resolution than the other stars in Figure 3. An absorption 
line of He I at A2.058 /xm is present for all three stars, although the detection is marginal 
for X174711.4. Weak N ill emission appears at A2.247 //m in the low-resolution spectra of 
X174711.4 and X174656.3, while this line appears as a strong doublet feature in the higher 
resolution spectrum of X174617.0. Weak N v emission is also observed in all three stars, 
although its detection in the higher resolution spectrum of X174617.0 is the most secure. 
X174617.0 also exhibits a weak He ll line (A2.189 iim) in absorption. 
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These features are consistent with 04-6If''" and core-H burning WNL stars of the 
hydrogen- rich variety (WNh, Smith et al. 1996). OIf+ and WN8-9h stars with relatively 
weak lines can be very difficult to distinguish from one another using near-infrared spec- 
troscopy alone, as these two spectral types occupy opposite ends of a continuous morpholog- 
ical sequence (Bohannan & Crowther 1999). WN8-9h stars are usually distinguished by He 
II (A2.189 /xm) in emission or as a P Cygni profile. However, several Olf"'" have been shown 
to exhibit this line in emission as well, while several WN8-9h stars are known to exhibit this 
fine in absorption (Conti et al. 1995; Bohannan & Crowther 1999). An additional classifi- 
cation criterion, more relevant for the K band, is provided by the Br7 emission line, which 
is usually stronger than the A2. 112-2. 115 //m complex for WN8-9h stars, relatively weak in 
Olf stars, and of comparable strength in Olf"'" stars; based upon this criterion, we classify 
X174656.3 and X174711.4 as WN8-9h stars, and maintain the original OGIf"*" classification 
from MMM07 for X174617.0. 

WN8-9h stars are also well represented in the Arches cluster. A model spectrum of the 
WN8-9h star F7 from Martins et al. (2008) is included in Figure 3 for comparison with the 
other WN8-9h stars. The basic parameters of this model are Tg// = 32, 900 K, L = IO^'^Lq, 
M = 2.5 X 10"^ Mq yr-^ and v^o = 1300 km s'^ 

Figure 4 presents the iT-band spectra of X174712.2, X174713.0, and X174516.7. These 
stars are distinguishable from the WN8-9h stars of Figure 3 by their more prominent He 
II emission at A2.189 /im, indicative of slightly earlier spectral types in the range of WN7- 
8h (e. g., sec Martins et al. 2008). The higher signal-to-noisc spectrum of X174712.2 also 
exhibits emission features of N v near A2.10 /xm, and C IV emission near A2.078 /^.m. In 
conclusion, we classify X174712.2, X174713.0, and X174516.7 as WN7-8h stars. 



3.4. Hydrogen-Deficient WN Stars 

Figure 5 presents K-hand spectra of WN stars X174550.6, X174508.9, X174555.3, and 
X174522.6. These stars have hydrogen-deficient spectra containing very broad, flat-topped 
emission lines of He I and He ll, indicative of fast, extended stellar winds. The discovery of 
X174550.6 was first reported in Cotera et al. (1999), and it was first identified as an X-ray 
source in Muno et al. (2006). We present a new, higher resolution X-band spectrum of this 
star and examine its spectral type. 

For X174550.6 we measure EW(A2.189 /xm)/EW(A2.1661 /xm)=0.52 and EW(A2.189 
/im)/EW(A2.112 /im)=0.55, which is consistent with stars of subtype WN7, according to 
Figer, McLean & Najarro (1997); this classification updates the earher WN6 classification in 
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Cotera et al. (1999). By direct comparison, the spectrum of X174508.9 is very similar to that 
of X174550.6, also exhibiting EW ratios consistent with WN7 classification. Slight differences 
in the spectrum of X174508.9 include a stronger He ll emission at A2.189 and A2.314 /xm, 
and a deeper P Cygni absorption component from He I at A2.058 fim. The He I absorption 
component for X174508.9 also has larger blueshift than that of X174550.6, implying that the 
wind of X174508.9 may be faster and more extended than that of X174550.6. 

The counterpart to X174555.3 was first presented in MMM07, where it was classified as 
a WN6b star. However, a refined spectral analysis of this star demonstrates that the original 
telluric correction and continuum fitting resulted in an underestimate of the He I (A2.0587 
/xm) emission-line flux. The more reliable result presented here exhibits a stronger emission 
line at this wavelength, with a blueshifted absorption profile from the same He I transition 
superimposed on the emission. This absorption line forms in the outer parts of the stellar 
wind where He ll has recombined, and so provides an estimate of the wind velocity there. 
We measure the absorption line centroid at A=2.0425 /xm, which implies a wind velocity of 
2400 ± 100 km s~^, typical for WN5-6 stars (Crowther & Smith 1996). Furthermore, the 
broad emission feature near A2.11 /xm appears to extend further into the blue end of the 
line, when compared with the same feature in the WN7 stars X174550.6 and X174508.9. We 
attribute this to N v emission, which is a characteristic of WN4-5 stars (Crowther et al. 
2006a). Therefore we relax our previous classification of WN6b to WN5-6b for X174555.3. 
Finally, we classify X174522.6 as a WN5-6b star as well, owing to the strong similarity 
between its spectrum and that of X174555.3. 



3.5. Weak-lined WC9 Stars 

Figure 6 presents the K-hand spectra of X174519.1, X174617.1 and X174645.2. The 
former two stars exhibit broad emission lines from the complex of He I, N ill, C ill and O 
III (A2. 112-2. 115 /xm). He I and He ll (A2.058 /xm), C ill and C IV (near A2.07-2.08 /xm), 
and weak He ll emission near A2.189 /xm. These spectral characteristics are consistent with 
those of late carbon-type (WC) Wolf-Rayet stars, specifically WC9 stars, whose values of 
EW(A2.08 /im)/EW(A2.112 /xm) are below 1 (e. g., see the spectrum of WR 112 in Figer, 
McLean & Najarro 1997). However, the peaks of the strongest emission lines are only 
~ 10% of the continuum fiux. WC9 stars typically exhibit much stronger lines. However, 
WC stars associated with hot thermal dust emission (WC9d stars) typically exhibit relatively 
weak emission lines in the A'-band, much like those of X174519.1 and X174617.7 (e. g., see 
Crowther et al. 2006a, their Figure 11). This is presumed to be the result of a bright thermal 
dust continuum which competes with the line emission of stellar wind. 
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The lower spectrum in Figure 6 is for XI 74645. 2, and appears completely featureless 
within the noise level. Based on this star's infrared SED and hard X-ray emission, it was 
identified by Hyodo et al. (2008) a dust-enshrouded, late-type WC binary, also known as a 
DWCL star (Williams et al. 1987), the same spectral type as the Quintuplet proper members 
(e. g., see Figer, McLean & Morris 1999). Figure 6 provides the first near-infrared spectrum 
of this star, and it is consistent with the DWCL spectral type based on its featureless 
character. Although the weakness or absence of emission lines in the near-infrared spectra 
of late-type WC stars is mainly attributable to dilution by continuum emission from hot 
dust, line dilution may also result from the stellar continuum emission from a very luminous 
companion. Later, in Section 4, we will present additional evidence for the presence of hot 
dust associated with X174519.1 and X174617.7. 



3.6. CXOGC J174516.1-290315 

X174516.1 is one of the brightest emission-hne stars known in the Galactic center (Kg — 
7.84 mag). This source was first identified as a young stellar object by FeUi et al. (2002), 
and it was subsequently characterized as an Ofpe or LBV emission-line star by Muno et 
al. (2006), who also showed this star to have a 2.3 mJy radio counterpart at uSA GHz. 
Stars of these spectral types are hydrogen rich, yet represent the beginning of post main- 
sequence evolution for the most massive stars (Morris et al. 1996). Figure 7 shows our 
i^T-band spectrum of X174516.1 at medium resolution. Strong emission lines of Br7 and He 
I are accompanied by weak low-ionization emission lines of Mg ll. The A2.058 /xm line is the 
strongest emission feature and appears as a P-Cygni profile, while the He I line at A2.112 /xm, 
also exhibiting a P-Cygni profile, is much weaker and may be blended with weak emission 
from N III. There also appears to be a weak P-Cygni feature of He ll near A2.19 //m. Based 
on the above characteristics, X174516.1 appears to be a member of Ofpe/WN9 class (e. g., 
see Morris et al. 1996), also resembling several stars that lie near the central parsec of the 
Galactic center, namely IRS16C, 1RS16NW, and IRS33E, the spectra of which were modeled 
in Martins et al. (2007). Owing to the strong similarity of our spectra to these stars, we 
classify X174516.1 as an Ofpe/WN9 star. 

From the P-Cygni profile of He I A2.058 /xm, we may infer the velocity at the radius 
in the wind where He I becomes an absorption line. The trough in the A2.0587 fim feature 
is centered at 2.0519 fim and the emission peak is centered at A2.0541 jj,m, so the implied 
velocity in the absorption zone is ~ 320 km s^^. 
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4. Stelleir Properties 

4.1. Bolometric Luminosity 

The near-infrared photometry of each massive star is presented in Table 2, which also 
includes mid-infrared photometry from the Spitzer Infrared Array Camera (IRAC) point- 
source catalog of the Galactic Center (Ramirez et al. 2008). To derive the absolute photom- 
etry, we must first determine the effects of interstellar extinction. Nishiyama et al. (2006) 
derived a near-infrared extinction relation for stars of the Galactic center using photometry 
from the SIRIUS survey; they report the following relations: A^^ = 1.44 ± 0.01Eh-Ks ^-iid 
Ak^ = 0.494 ± 0.006^j_x,, where E{H - K,) = {H - Ks)ohs - {H - Ks)o. For the stars 
presented in this work, the observed JHKg photometry was generally taken from the SlRlUS 
catalog, except for very bright stars that were saturated {Kg < 8 mag), in which case we 
used 2MASS photometry. The intrinsic colors, {H — Ks)q and (J — Ks)q, were taken from 
the hterature using comparison stars of similar spectral types. For the WN stars, we used 
the intrinsic colors from Crowther et al. (2006a, their Table Al). For the O supcrgiants we 
used the synthetic photometry for O supcrgiants from Martins & Plez (2006). Using the 
adopted intrinsic colors, we were able derive two values of Ak^ using the observed J — Ks 
and H — Kg colors, which were then averaged to produce a final value for the Ks-band 
extinction. The absolute i^-band magnitude of each star was derived assuming a distance 
of 8 kpc to the Galactic center. To calculate the total luminosity of each star, bolometric 
corrections of comparison stars having similar spectral types were also taken from the same 
literature sources referenced above. In the interest of a thorough comparison, we not only 
performed these calculations for new WN and O X-ray sources, but included previously 
known sources as well. The results arc presented in Table 3. The most luminous stars in 
the sample arc the WN8-9h stars X174656.3 and A7 (Figer et al. 2002) of the Arches clus- 
ter (CXOGC J174550.4-284919), with Lboi = 2.9 x 10^ Lgand 3.9 x 10^ Lq, respectively; 
the least luminous star is the broad-lined WN7 star X174508.9, which has a luminosity of 
Lboi = 2.5 X 10^ Lq. In contrast to WNh stars which are still burning hydrogen, the WN7 
stars, such as X174508.9, are in an evolved state of core-He burning. Such stars have already 
blown off their hydrogen envelope, and as a result, have a significantly smaller radius and 
lower bolometric luminosity than the less-evolved WNh sources. 

The situation is more difficult in the case of the WC stars; we were unable to derive 
an intrinsic color for these stars, owing to their significant near-infrared excess. Instead, the 
photometry of X174617.7 was extinction-corrected by adopting the average extinction value 
for the Quintuplet cluster [Ak — 2.7 mag; Figer et al. 1999), since X174617.7 appears to 
be an outer member of the Quintuplet. For X174519.1, wc adopted the extinction value 
measured for a nearby B2Ia star, S174523.1 from MMM07, which has Ak — 3.2 mag. 
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Although these extinction approximations will be useful in the analysis of the WC star 
SEDs, we did not compute bolometric luminosities for them, and they are not included in 
Table 3. 



4.2. Infrared Excess 

To determine whether the massive stellar X-ray sources are distinguishable from Galactic 
center field stars in the infrared, we constructed color-color diagrams. Figure 8 illustrates 
that almost every star exhibits a significant infrared excess, apparent as an offset from the 
main reddening locus that is traced by the roughly linear scatter distribution of field stars. 
In H — Kg versus J — H color space, the WN stars and most of the O supergiants occupy the 
same space, while the WC9 and DWCL stars appear more separable as a group, owing to 
their significantly larger i^^-band excess. With the inclusion of A8 /xm photometry, the WN, 
O, WC9, and DWCL stars are even easier to distinguish as separate groups in Kg — [8.0] 
versus J — Kg color space. The WN stars and one O supergiant appear to cluster together, 
while the WC stars, again, exhibit a relatively large excess compared with the WN stars, 
especially the DWCL stars qF211 and qF231 of the Quintuplet cluster. Three O supergiants 
exhibit an extremely large excess at A8 ;um, which likely originates from unresolved, warm 
gas and dust surrounding these young stars (evidence in support of this hypothesis provided 
by mid-infrared field images is presented in Section 6.2). 

There are several ways in which Wolf-Rayet and O stars can exhibit an infrared excess. 
Within the ionized winds of hot supergiants, scattering of electrons off of H"*" and He"'" 
ions produces free-free emission, which contributes a significant amount of continuum flux 
in the infrared, becoming more significant with increasing wavelength as the underlying 
photospheric continuum becomes fainter with respect to the free-free component (Wright & 
Barlow 1975). Thermal emission from hot circumstellar dust (T = 600-1300 K) is another 
source of infrared excess that is often associated with WC stars (Cohen et al. 1975, Williams 
et al. 1987). This form of excess emission can exhibit significantly higher flux density than 
free-free emission. Finally, a large infrared excess may also be detected from massive stars 
surrounded by unresolved clouds of gas and dust local to the star. 

To further investigate the nature of the infrared excesses observed in our sample, we 
generated Al-10 //m SEDs. For the photometry of the WN and O stars we applied extinction 

corrections to the photometry in each band according to the extinction ratios in Indebetouw 
et al. (2005), using the derived values of in Table 3 and the adopted values for the 
WC stars. To convert the mid-infrared data to flux density units, we applied zero points 
of 280.9, 179.7, 115.0, and 64.13 Jy for IRAC channels 1, 2, 3, and 4, respectively (IRAC 
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Data Handbook-Table 5.1). In the absence of confusion-induced photometric errors, the 
near-infrared SIRIUS photometry matches 2MASS very closely; the mean and rms variance 
of the magnitude differences are 0.009 and 0.015 mag in the J band, respectively (Nishiyama 
et al. 2008). Therefore, we obtained near-infrared fluxes for SlRlUS using the 2MASS zero 
points of 1594.0, 1024.0, and 666.7 Jy for the J, H , and Kg bands, respectively. Below, we 
examine the SEDs of several WN and WC stars to investigate the nature of the infrared 
excess emission in more detail. 



4.3. Free-Free Emission from WN/O Stairs 

To test the hypothesis that free-free emission is responsible for the infrared excess ob- 
served toward the massive stars in our sample, we employed a modifled version of the free- free 
emission model of Cohen et al. (1975). In this model the optically-thin, free-free emission 
spectrum for an unresolved source may be written as 



F,(Jy) = 2.3 X 10-^^ :g^!£^e-c/ATe 

where Z is the mean electric charge per ion, and Ni are the respective electron and ion 
densities in units of cm~^, Tg is the electron temperature in Kelvins, Rg is the size of the 
emitting region in cm, d is the distance to the source in cm, A is the wavelength in /xm, 
and C= 14400 //m K. Free- free emission from evolved massive-star winds typically becomes 
optically thick in the mid-infrared at some wavelength, Aj. If this turnover wavelength is 
observed, one can also solve for the emission measure: 



N N R = 7.3 X 10=^^ y/?; 5 

iVeJv,/t, AfZ2(l - e-^A*^«) ■ ^ ' 

By combining Equations (1) and (2), we constructed model SEDs for several stars: 
X174516.1(Ofpe/WN9), X174555.3 (WN5-6b), X174656.3 (WN8-9h), and X174712.2 (WN7- 
8h), all of which have an observable turnover wavelength in their mid-infrared SEDs. For 
each star we determined the flux density attributable to free-free emission by subtracting 
off a blackbody continuum. The continuum flux was calculated using the luminosity values 
in Table 3 and the Stefan-Boltzmann relation L = AnR'^aT'^, where L is the bolometric 
luminosity, R is the stellar radius, and T is the effective temperature. Unfortunately, there 
is a degeneracy between temperature and radius on the Rayleigh-Jeans portion of the SED, 
which we are Umited to. So, initially we guessed values of T and R that are appropriate for 
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the spectral type of the star being analyzed. The values were then adjusted until we ob- 
tained a satisfactory match to the extinction-corrected near-infrared photometry. Of course, 
a slight underestimate of T could be compensated by an overestimate of R. This uncer- 
tainty, however, is unimportant for our purposes because the free-free emission spectrum is 
not affected by it. Once the resulting Rayleigh- Jeans component was subtracted from the 
photometry, the remaining free-free flux, described by Equations (1) and (2), was fit using 
a least squares method. By performing a least squares fit to the excess spectrum, we were 
ultimately able to derive Te, TVe, At, and Rg. For the models we assume A^e — and Z — 1 
for the WNh and Ofpe/WN9 stars, and Z = 2 and Ne = 2Ni for the hydrogen-deficient 
WN5-6b star. 

The resulting model parameters are listed in Table 4, and are illustrated graphically in 
Figure 9. The model curves match the data reasonably well, and the derived values for Tg, 
Ne, and Rg are similar to those of the WN stars modeled in Cohen et al. (1975). Rehable 
errors in our parameter fits are difficult to ascertain; uncertainties in the extinction law and 
the effect of strong emission lines on the photometry are likely to play a role. Nonetheless, 
a very precise free-free emission model of the stars in our sample is not the focus here. The 
goal is to unambiguously demonstrate that the infrared excess emission from the WN stars 
in Figure 8 is well described by a simple free-free emission model having physical parameters 
that are typical for such stars. 



4.4. Hot Dust Emission from WC Stars 

The infrared excess emission from the WC stars in our sample (X174519.1 and 174617.7) 
has a significantly higher intensity than that of the WN stars, which is apparent in Figure 8. 
Figure 10 shows the extinction-corrected SEDs of the new WC9 discoveries X174617.7 and 
X174519.1, accompanied by the SEDs of the Quintuplet proper member X174614.6 (qF231 
fi:om Figer, McLean & Morris 1999), a putatively single WC9 star HDM13 (Hadfield et al. 
2007), and the WC9d star WR 59 for comparison. For the photometry of the comparison 
stars HDM13 and WR 59, we applied extinction corrections using the Ak or Ay values from 
Hadfield et al. (2007) and van der Hucht (2001), respectively. 

For the WC9 stars X174617.7 and X174519.1, we were unable to use our free-free emis- 
sion model to simultaneously reproduce the near-infrared and mid-infrared fiux densities. 
No sensible combination of the stellar radii, effective temperature, and wind parameters (Tg 
and Rg) resulted in a reproduction of the rise in fiux density between the H and Kg bands. 
Fitting the bright mid-infrared portion of the WC9 SEDs without overestimating the near- 
infrared portion would require insensibly low values of Tg. So, even though WC stars are 
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indeed important sources of infrared free-free emission, a stronger additional emission com- 
ponent is required to explain the structure of their infrared SEDs. Thermal emission from 
hot dust has been observed from many WC9 stars and binaries, and we suspect this is the 
case for X174617.7 and X174519.1. This hypothesis is strengthened by the following compar- 
ison with the other dusty WC stars in Figure 10. qF231 exhibits the largest infrared excess. 
This star is an X-ray source (X174614.6), and has been classified as a DWCL star (Figer et 
al. 1999; Moneti et al. 2001). Such stars exhibit nearly featureless near-infrared spectra, as 
demonstrated in Section 3.5 (Figure 6) with the spectrum of the DWCL star X174645.2 (the 
Spitzer /IRAC photometry of X174645.2 is saturated, so its SED was not included in Figure 
10). By comparison, the infrared excess from the single WC9 star HDM13 is the smallest of 
the stars in Figure 10, exhibiting characteristics very similar to the WN stars in Figure 9. 
Thus, the excess from HDM13 can perhaps be attributed to free-free emission alone. The 
magnitude of the infrared excess in the SEDs of X174519.1 and X174617.7 appears to lie 
in-between the extremes defined by the SEDs of qF231 and HDM13, but bear the closest 
resemblance to the WC9d star WR 59. Compared with qF231, the smaller infrared excess of 
X174519.1 and X174617.7 provides a natural explanation for why the emission lines in their 
spectra are only partially diluted, instead of being completely diluted as in the case of the 
DWCL spectra. Still, the mid-infrared SEDs of these sources exhibit an excess component 
that rises significantly farther above the RJ continuum than the free-free components of the 
single WC9 star HDM13 and the WN stars in the previous section. Furthermore, the excess 
emission peaks between A5 and AlO fim and is detectable at wavelengths as short as the 
J-band wavelength. Thus, these characteristics can be attributed to the presence of hot 
dust with temperatures near ~ 1000 K. Realistically, the collective excess of these objects 
is probably due to multiple hot dust components lying at different radii from the central 
star(s), stellar continuum photons, and free-free emission. A more detailed modeling of the 
complex excess emission for these stars is beyond the scope of this paper. We simply intend 
to demonstrate that hot dust is responsible for the bulk of the mid-infrared excess detected 
from X174519.1 and X174617.7, and provides a simple explanation for the continuum diluted 
emission lines observed in their i^-band spectra. 

WC stars have long been known to be producers of dust if they exist in a binary 
system. In fact, it has been suggested that all dusty WC stars (WCd) are members of 
binaries (Williams et al. 2005), in which the dust is produced by the collision of carbon- 
rich and hydrogen-rich winds from the WC star and an OB companion, respectively. Wind 
collision also naturally explains the detection of these sources at hard X-ray energies (see 
Section 6.1.1). For several known WCd stars, the binary hypothesis has been strengthened 
beyond any reasonable doubt by high-resolution imaging experiments, which have revealed 
rotating Archimedian spirals of dust surrounding them (e. g., seeTuthill et al. 2006, 2008). 
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Similar structures may be observable, modulo unfortunate projection effects, in association 
with X174519.1, X174617.7, and X174645.2. 

5. X-ray Emission 

Table 5 lists the X-ray photometry for all 30 confirmed massive stellar X-ray sources in 
the Galactic center region, including new and previously identified objects. The data for each 
star were compiled from the catalog of Muno et al. (2009). Table 5 includes X-ray source 

astromctry, total on-source integration time, photon number counts in the hard and soft 
bands, broadband photon flux, hardness ratios for soft and hard energy bands, the average 
energy of all detected photons, the derived X-ray luminosity, and the spectral types of their 
stellar counterparts. Most of the sources have photon fluxes ranging from Fx = 2 x 10^'' to 
Fx = 7 X 10^^ cm~^ s^^. An approximate energy flux can be obtained by multiplying Fx by 
the average photon energy (F), which yields values ranging from 1.3 x 10~^^ to 5.7 x 10~^^ 
erg cm~^ s~^. 

In order to constrain the nature of the X-ray emission, it is useful to examine the 
hardness ratio, defined as HR=(/i — s)/{h + s), where h and s are the fiuxes in the hard 
and soft energy bands, respectively. The soft color, defined as HRO, is defined by letting h 
and s be the fiuxes in the respective 0.5-2.0 keV and 2.0-3.3 keV energy bands. In Paper I, 
we used HRO to select X-ray sources that are hkely to be located near the Galactic center, 
since the softest X-ray photons are highly absorbed by the ISM. For this reason, many of 
the sources in Table 5 are not detected in the softest 0.5-2.0 keV energy band, which results 
in them having an HRO value of 1. The hard color, HR2, is defined by letting h and s be 
the fluxes in the 4.7-8.0 keV and 3.3-4.7 keV energy bands, respectively. Compared with 
HRO, the more energetic photons used to calculate HR2 suffer relatively little absorption 
from intervening gas and dust, so they are more useful for constraining the hard components 
of the X-ray SEDs. 

Figure 11 is a scatter plot of HR2 versus broadband (0.2-8.0 keV) photon flux for 
absorbed X-ray sources that are likely to be located near the Galactic center (i.e., those 
sources whose soft X-ray colors imply a hydrogen absorption column of > 4 x 10^^ cm~^). 
The majority of the X-ray sources that comprise the "cloud" of small data points in Figure 
11 are hkely to be CVs, a hypothesis that is consistent with both the X-ray properties 
of the sample (Muno et al. 2003, 2009) and the rarity of near-infrared counterparts with 
Ks < 15.6 mag, discussed in Paper I. By comparison to the fleld population of presumed 
CVs, the conflrmed massive stars appear to be systematically softer and brighter in X-rays. 
To estimate the plasma temperatures of thermal sources, or the power-law indices of non- 



- 18 - 



thermal sources, we constructed models of HR2 versus luminosity for optically-thin thermal 
plasmas with temperatures of kT = 0.9, 2.2, 3.4, and 6.8 keV, and power-law emission 
models for non-thermal sources with photon indices of — 1 < F < 3. The models were 
constructed for sources lying at a distance of 8 kpc that are absorbed by a hydrogen column 
of iVn = 6.0 X 10^^ cm~^. This value of A^h is generally adopted as a reasonable average 
value for X-ray sources near the Galactic center (Muno et al. 2003), and is also consistent 
with the average extinction of the stars listed in Table 3, according to the relation between 
A'h and Ay given in Predehl & Schmitt (1995). All models were computed using the XSPEcE 
program. The HR2 values of the massive stars imply the presence of plasma components 
with temperatures of kT > 1 keV (in most cases kT > 2 keV), if the emission is thermal, 
and luminosities in the range ~ lO'^^-lO^^ erg s~^ (0.5-8.0 keV). The X-ray sources with 
the highest fluxes and hardest colors are the DWCL star X174645.2 (Hyodo et al. 2008) and 
several of the Arches cluster members, all of which appear to have relatively high luminosities 
near ~ 10^^ erg s~^. The softest sources, undetected in the 4.7-8.0 keV band, are the 
O9I-B0I star X174537.3 and the Quintuplet members qF231 (DWCL star) and qF242 (OBI). 
The hardest massive stellar X-ray source in the Galactic center is the WN7 star X174508.9, 
with HR2=0.3, and an average photon energy of 5.3 keV. We note that a constant value of 
A'h is a gross simplification; more so is the assumption of a single temperature component 
for the X-ray plasma. Below, we address the effect that these simplifications have on the 
derived properties of the X-ray-emitting plasma. 

The fiducial conversion from photon flux to unabsorbed 0.5-8.0 keV luminosity is such 
that 10^^ erg s~^ equals 6 x 10~^ cm~^ s~^, assuming F = 1.5 or kT = 7 keV, which are 
average values for the entire X-ray population. Assuming kT = 2 keV instead, which is 
more applicable to our sample of massive stellar X-ray sources, the conversion is such that 
Iq3a gj,g g-i gq^aig 4 x 10~^ cm~^ s~^. However, the uncertainty in these conversion factors 
is substantial. Indeed, interstellar absorption has extinguished most of the X-ray flux at 
energies below 1 keV, which is where massive stars emit the bulk of their X-rays (e.g., see 
Oskinova 2005, Skinner et al. 2002a), so we could be significantly underestimating the total 
intrinsic X-ray luminosities in the 0.5-8.0 keV range. As a particularly fitting example of 
the errors imposed by assuming a single temperature plasma, we refer to the X-ray analysis 
of the source X174536.1 by Mikles et al. (2006). Owing to relatively bright soft and hard 
X-ray fluxes, the authors were able to fit a detailed thermal plasma model to the X-ray 
spectrum of this source. Approximating with a single temperature plasma, as we do for 
our sample as a whole, the authors derived kT ^ 3.5 keV and Lx = 10^^'^ erg s~^ (0.5-8.0 
keV) for X174536.1, which is reasonably close to the plasma temperature implied by the 



^http: / /heasarc. gsfc.nasa.gov/docs/xanadu/xspec/ 
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HR2 value of this source, and to our derived luminosity of 10^^-^ erg s"^. However, the 
authors obtained a more accurate spectral fit using a two-temperature model, having a very 
soft component kT ^ 0.7 keV and an additional hard component with kT ^ 4.6 keV, each 
with respective luminosities of = 10^^-^ erg s^^ (0.5-2.0 keV) and 10^^-^ erg s^^ (2.0- 
8.0 keV). Thus, the resulting X-ray luminosity difference between the 1-T and 2-T plasma 
models is a sobering pa 1.4 dex. Clearly, in assigning a single temperature plasma we may be 
neglecting a significant amount of soft X-ray luminosity, if it lies beneath our fiux threshold 
due to absorption. Unfortunately, only two of the new detections reported in this work 
were detected at energies below 2 keV, the brightest of which has only 13 photon counts in 
that range. Without a significant fiux of soft X-rays to subject to a model fit, assuming a 
single temperature plasma is the best we can currently do. So, we assumed that 2 keV is 
a reasonable average intrinsic plasma energy and based our luminosity calculations on this 
assumption. The resulting X-ray luminosities are included in Column 12 of Table 5. As a 
means of estimating the errors, we compared our values to the results of other authors who 
were able to derive luminosities for relatively bright. Galactic center X-ray sources based 
on a modeling of their X-ray spectra, which enable them to solve for kT directly. We did 
not use the above-mentioned X174536.1 in this exercise because it is a variable X-ray source 
(see Section 5.1). The comparison is facihtated by Table 6, which includes the three bright 
WN8~9h X-ray sources from the Arches cluster (Wang, Dong & Lang 2006), the DWCL star 
X174645.2 (Hyodo et al. 2008), the WN7 star X174550.6, and the Of star X174528.2 (Muno 
et al. 2006). The average and standard deviation of the difference in the X-ray luminosity 
between our results and those of the other authors is 0.25 (0.24) dex. Thus, we assign 3a 
errors of 0.75 dex for our derived luminosities. The results should still be met with caution, 
however, since the comparison measurements from the other authors listed in Table 6 also 
assumed a single temperature plasmas. Finally, using our derived X-ray luminosities and 
the stellar bolometric luminosities in Table 3, we recovered a relation between these values, 
such that Lx/Lb„i = iQ-e-Sio.e pigure 12 illustrates this result. 

In conclusion, the main points we intend to convey with X-ray data are that (1) the 
X-ray photometry of the massive stars as a group is systematically softer and brighter than 
the field population presumed to be dominated by CVs, (2) each source contains a hard 
X-ray component, requiring either a thermal plasma with kT > 1 keV (in most cases kT > 2 
keV), or a power-law source with —1 < F < 3, (3) the majority of sources have 0.5-8.0 keV 
luminosities between ~ 10^^ and 10^^ erg s~^, while several others appear to have higher 
luminosities closer to 10^^ erg s^^, and (4) the Lx/I/boi ratio for the Galactic center sources 
has an approximately constant value of io~^-^='=°-^. 
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5.1. X-ray Variability 

Variability of the Galactic center X-ray population was addressed in Muno et al. (2009). 
Variable sources were grouped according to three types of variability: flux variations between 
individual observations, separated by day- to year-long time scales (long-term variability), 
variations that occurred within individual observations (short-term variabihty), and periodic 
variations that occurred within an individual observation. Four of the massive stellar X-ray 
sources near the Galactic center exhibit long-term variability. None of them exhibits short- 
term or periodic variability. The variable sources include the WN7-8h stars X174536.1 and 
X174712.2, the Ofpe/WN9 star X174516.1, and the OBI Quintuplet member X174614.5. We 
reproduce the variability data for these sources from Muno et al. (2008, their Table 6) in our 
Table 7. The table includes the source record numbers from the original X-ray catalog, the 
observation ID numbers in which the largest and smallest fluxes were observed, the values 
of the largest and smallest fluxes, and the ratios of those fluxes. The most variable source 
among the massive stars is the WN7-8h star X174712.2, which varied by a factor of 3.6 
between 2000 March 29 and 2001 July 16. 

A variety of phenomena may give rise to X-ray variability. GoUiding-wind binaries 
with eccentric orbits will increase in X-ray luminosity near periastron. Eclipses of a wind 
collision zone may also result in variable fluxes if the orbit is close to edge-on. However, 
colliding- wind binaries may not exhibit detectable X-ray variations, as in the case of WR 25 
(Pollock & Gorcoran 2005), which remained stable for a time span of ten years (Raassen et al. 
2003). For accreting compact objects in HMXBs, eclipses are also a possible soiuxe of X-ray 
variability, as is an intrinsic variation of accretion rate. Finally, X-rays from hot spots within 
an equatorial disk surrounding a massive star (see Section 6.1) may be modulated with the 
stellar rotation period. The applicability of these models in explaining the characteristics of 
the X-ray emission of our sample is discussed in the following section. 

6. Discussion 

6.1. Physical Origin of the X-ray Emission 

The X-ray photometry of the majority of massive stellar X-ray sources near the Galactic 
center indicates that the emission is from thermal plasmas having temperatures of kT > 1-2 
keV, or non-thermal sources with power-law indices of — 1 < F < 3, and X-ray luminosities in 
the range Lx = 10^^-10^^ erg s~^ (0.5—8.0 keV). Hard- X-ray emission at these luminosities 
is not a ubiquitous feature of single O supergiants or Wolf-Rayet stars; these features are 
typically regarded as an indication of binarity, since it is well established that hard X-rays 
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can be generated within the shocked interface of the binary components' opposing supersonic 
winds, or by accretion of stellar wind material onto a compact companion. However, there 
are extraordinary mechanisms by which some single massive stars can intrinsically generate 
hard X-rays as well; so, we entertained several interpretations of the hard X-ray emission 
from our sample. 

6.1.1. Intrinsic X-ray Emission from Massive Stars 

The intrinsic X-ray emission from single massive stars throughout the Galaxy is typically 
thermal and soft, having kT ~ 0.6 keV (e.g., see Oskinova 2005 and references therein). 
The X-ray luminosities of O stars follow an observed trend in which they scale with stellar 
bolomctric luminosity as Lx/-f/boi ~ 10~^ (Long & White 1980; Bcrghocfcr ct al. 1997; 
Sana et al. 2006; Broos et al. 2007, Albacete Colombo et al. 2007). This trend has been 
interpreted as a manifestation of the physical link between stellar luminosity, which drives 
supersonic winds, and the X-ray-emitting shocks that result from those winds via line-driven 
instabilities (Lucy & White 1980; Owocki et al. 1988; Feldmeier, Puis, & Pauldrach 1997). 
In Section 5, we derived the average and standard deviation of the bolometric to X-ray 
luminosity ratio for our sample, and obtained Ly^/L^^^ = iQ~^-^^o-^^ which is consistent with 
the known trend for single massive stars. However, the fact that the Lx/i^boi relation is 
obeyed does not necessarily imply that the X-rays arc solely generated via radiative shocks 
within a single stellar wind; colliding- wind binaries have also been shown to obey the Lx/I>boi 
relation as well, although many colliding- wind binaries have Lx/I>boi closer to 10~^ (e.g., see 
Oskinova 2005). Still, the emission of hard X-rays is uncommon property of single massive 
stars. There is, however, an extraordinary mechanism by which massive stars can genrerate 
hard X-rays intrinsically. The extraordinary 05.5V star 9^ Orionis C, for example, has an 
average X-ray luminosity of Lx ~ 10^^ erg s~^ (0.5-10.0 keV), and a hard spectral component 
from thermal plasma with kT ^ 3 keV (Gagne et al. 2005). The hard X-ray emission is 
modulated with the star's 15 day rotation period; so, the hard X-rays have been interpreted 
as a result of the star's strong, fossil magnetic field {B ^ 1000 G, Donati et al. 2002), which 
channels and confines the outflowing stellar wind into an equatorial collision/cooling disk 
(the magnetic wind-confinement model) . The process results in strong shocks that generate 
plasma temperatures in excess of 10^ K, and hence, hard X-rays. The pecuhar Of ?p star HD 
191612 is another example, having a measured magnetic field of 1500 G (Donati et al. 2006) 
and an X-ray spectrum containing thermal plasma components with kT ^ 1-2.5 keV and 
an X-ray luminosity of 7-9 x 10^^ erg s ~^ (Naze et al. 2007). HD 191612 exhibits periodic 
spectroscopic variations (a time-variable equivalent width of the Ha line) with a period of 
538 days. Walborn et al. (2004) discussed the possibility that the variability is induced 
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by the periastron passage of a binary companion, while Donati et al. (2006) attributed 
the variations to the star's rotational period and, hence, its strong magnetic field. Since 9^ 
Orionis C and HD 191612 both exhibit periodic changes in their spectra, we should consider 
the potential to observe the same phenomena in the spectra of the O supergiants of our 
sample, if they represent similar stellar phenomena. 

The magnetic confinement model has been discussed in attempts to explain the detec- 
tion of hard X-rays from the nitrogen-rich Wolf-Rayct stars as well, in particular, for WRllO, 
EZ CMa, and the oxygen-rich Wolf-Rayet star WR142, all of which have stellar and X-ray 
properties that are very similar to our sample and exhibit no evidence for a companion (Skin- 
ner et al. 2002a, 2002b; Oskinova et al. 2009). However, the magnetic wind-confinement 
model faces difficulties when applied to Wolf-Rayet stars; the wind momenta of Wolf-Rayet 
stars are over an order-of- magnitude stronger than that of O stars such as 9^ Orionis C and 
HD 1919612, and are more difficult to channel effectively. According to Babel & Montmerle 
(1997), a stellar wind is confineable if the magnetic energy density exceeds the wind kinetic 
energy density {B'^/4tt > pv"^). Assuming p = M/Anvr'^, and v — O.bVoo at r = IRgtar, this 
would imply that magnetic fields greater than ^5 kG in strength are necessary to confine the 
winds of the Wolf-Rayet stars and O supergiants that dominate our sample. Currently, there 
is no evidence for such strong fields in Wolf-Rayet stars. Strong fields have been considered in 
the interpretation of hard X-ray emission from the WO star WR 142 (Oskinova et al. 2009); 
however, this source is intrinsically ~10 times fainter in X-rays than the faintest source in 
our sample. Thus, we do not currently regard the magnetic wind-confinement scenario as a 
highly viable model for Galactic center Wolf-Rayet stars, although it may explain the X-ray 
emission from some of the O stars in our sample. Still, we can not completely rule out the 
possibility that Wolf-Rayct stars can generate hard X-rays intrinsically by some other exotic 
mechanism, such as inverse Compton scattering (e.g., see Chen & White 1991). However, 
since there arc so many Wolf-Rayet and O stars known in the Galactic center that have not 
been detected in X-rays (e.g., in the Arches and Quintuplet clusters), intrinsic hard X-ray 
generation from our sample, which represents the minority fraction of known massive stars 
in the Galactic center, would imply that our particular stars are special, having physical 
properties not ubiquitous to all Wolf-Rayet stars and O supergiants in the region, which 
would be puzzling. 

6.1.2. Colliding- Wind Binaries 

The collision of opposing, supersonic winds in massive binaries creates a hot shock 
{T > 10^ K) that is a potential source of hard thermal X-rays with kT > 1 keV. It is 



-23- 



reasonable to adopt this model in interpreting the hard X-ray emission from our sample, 
since many known colliding-wind binaries exhibit hard X-rays with luminosities that are 
comparable to our Galactic center sample, such as WR 147 (WN8-I-0B; Skinner et al. 2006), 
WR 25 (WN6ha+04f ; Raassen et al. 2003; Albacete Colombo ct al. 2008), 7^ Velorum 
(WC8+07.5; Skinner et al. 2001), V444 Cyg (WN5+06; Maeda et al. 1999), Cyg 0B2 8 
(06If+05.5III(f); Albacete Colombo et al. 2007), and Plaskett's star (061+07.51; Linder 
et al. 2006), all of which have hard spectral components with kT ?a 2-4 keV and Lx = 10^^- 
lO^'' erg s-^ With regard to the WC stars of our sample (X174519.1, X174617.7, and 
X174645.2), the case for colliding winds is particularly strong, since most single WC stars 
are not detected in X-rays at all (Skinner et al. 2006), and hot dust is present (see Section 
4.3). It has been proposed that the production of dust from carbon-rich WC winds requires 
collision and mixing with a hydrogen-rich companion star's wind (Crowther 2003). While 
the explanation for hard X-ray emission in the WN/0 stars of our sample is not as clear, 
the colliding-wind hypothesis may be strengthened by the fact that binarity is common in 
Wolf-Rayet and O stars (Wallace 2007). 

If the colliding-wind binary hypothesis is correct, one might expect the signature of a 
companion star to be evident in the infrared spectra, presented in Section 3. However, the 
current infrared data are probably of insufficient spectral resolution and signal to noise to 
distinguish the contributions from individual stellar components of a massive binary. In 
the case where a Wolf-Rayet star or O supergiant has a fainter OB dwarf companion, the 
magnitudes may differ by a factor of ~3. With the implied flux ratio of Ril6 for the stellar 
components, H and He emission or absorption lines from the OB companion will be difficult to 
detect in a moderate-resolution spectrum. In the case where the binary consists of two similar 
Wolf-Rayet stars of comparable brightness, the broad emission lines of both components may 
blend together. In this case, the spectral lines of individual binary components might only be 
distinguishable with multi-epoch observations, where radial velocity shifts could betray their 
presence. Thus, the colliding-wind binary scenario is consistent with the current infrared 
spectra. 

Assuming our sample is mainly represented by binaries, it is worthwhile to speculate 
whether their global X-ray properties can be used to constrain the physical parameters of 
the binary system. For colliding-wind binaries, the emergent X-ray luminosity from the 
shocked wind-wind interface is a function of the wind velocities and mass-loss rates of the 
stellar components, and the binary separation. The exact functional form depends on the 
nature of the wind collision. For instance, if the wind momenta of two stars in a binary are 
comparable, then the shock zone will lie somewhere near the center of the space between 
them, where the opposing winds achieve momentum balance. The collision interface will 
form a bow shock structure that curves toward the star with the weaker wind. According 
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to Usov (1992; see their Equations (89) and (95)), the total X-ray luminosity in this case 
can be expressed as the sum of the contributions from two separate layers, one that lies 
on the inside (concave) surface of the bow shock, composed of material from the wind 
of the weaker star, and one that lies on the outside (convex) surface of the bow shock, 
composed primarily of material from the wind of the stronger star. Using mass-loss rates 
and terminal wind velocities that are typical for the Wolf-Rayet and O stars of our sample 
(M ~ IO-^-IO-^Mq yr-\ v^o = 1000-2000 km s-^ see Section 3), this particular model 
predicts X-ray luminosities of Lx ~ 10^^ — 10^^ ergs for binary separations oi D — 0.1-10 
AU. Given our uncertainty in X-ray luminosity, the brightest several sources in our sample 
could be binaries in which the components have comparable wind momenta. However, the 
majority of our sources have fainter X-ray luminosities in the range Lx ~ 10^^ — 10^^ erg s~^. 
The emergent X-ray luminosity in this particular colliding-wind model has a dependence of 
Lx oc D^^, so, this scenario could result in X-ray luminosities consistent with those of our 
sample only if the binary separations are very large {D > 10-500 AU). This would make the 
Galactic center sources analogous to the WN8h+0B binary WR 147, which has an estimated 
binary separation oi D 400 AU, and an X-ray luminosity and plasma temperature of 
Lx(0.5 - 7 keV) ^6x 10^^ erg s'^ and A;^ = 2 - 4 keV, respectively (Skinner et al. 2007). 
With periods of ~1000 years, binaries like these would not exhibit detectable radial velocity 
variations. However, at a distance of 8 kpc, they could have apparent separations of ~50 
mas, which could be resolvable with adaptive optics on 8-10 m class telescopes. 

Alternatively, Usov (1992) also considers the case where the ram pressure of the stronger 
stellar wind suppresses the weaker wind of the companion (see their Equation (81)), in which 
case the companion star is modeled as a hard sphere and the wind collision occurs very near 
or on its surface. For this particular model, which has a much stronger dependence on the 
binary separation (Z)^^), the expected X-ray emissivity of the shock zone yields lower values 
of X-ray luminosity in the range of ~ 10^^-10'^^ ergs s^^. Thus, if this particular model is 
applicable, then perhaps the Galactic center sources are binaries where the components have 
significantly different wind strengths. This model has been considered to explain the X-ray 
emission of the WN5-6 star WR 110 (Skinner et al. 2002a), which has Lx(0.5-7 keV) ~ 10=^^ 
erg s~^ and a hard spectral component with kT >3 keV. 

It is interesting that the most X-ray-luminous stars in our sample, with Lx ~ 10^^ 
erg s~^, are mainly comprised of objects within the Arches cluster. These particular stars 
could be binaries that have relatively small separations. Frequent stellar encounters within 
dense stellar clusters will tend to harden close binaries while 'ionizing' wide ones; so, it 
is reasonable to speculate that the relatively high X-ray luminosity of the Arches X-ray 
sources is the result of their relatively small orbital separations that were imposed by the 
dense cluster environment. If so, radial velocity variations will occur on shorter timescales 
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than for the rest of our sample and may exhibit a larger amplitude in the radial velocity 
variations, similar to that of the WN6h+0 binary NGC 3603-Al, which lies near the center 
of the starburst cluster NGC 3603 and exhibits radial velocity amplitudes of ~ 400 km s~^ 
(Schnurr et al. 2008). 



6.1.3. Supergiant High-Mass X-Ray Binaries 

We must also consider the possibihty that wind- accreting neutron stars or black holes 

are responsible for the X-ray emission from some of the massive stellar X-ray sources near 
the Galactic center, making them supergiant X-ray binaries (SGXBs). Although there is 
no reason to expect such objects to follow the observed trend of Lx/I>boi = 10~^, there are 
several quiescent SGXBs whose Lx/i^boi ratios that are consistent with this trend. So, the 
wind-accretion model still deserves consideration. 

According to Pfahl, Rappaport & Podsiadlowski (2002), the time-averaged 1-10 keV 
accretion luminosity of a wind-accreting neutron star is given by 
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e^)-^/^ (3) 



where 1] is the conversion efficiency of gravitational energy into radiation, M^, and are the 
mass-loss rate and wind velocity of the donor star, respectively, a is the orbital separation 
in AU, and e is the orbital eccentricity of the system. This formula is based on the Bondi- 
Hoyle-Lyttleton accretion model (Hoyle & Lyttleton 1941; Bondi & Hoyle 1944). For the 
massive stars of our sample, most of which have M ~ 10~^-10"^ Mq yr~^ and v ^ 2000 km 
s^^. Equation (3) implies that the emergent X-ray luminosity from a wind-accreting SGXB 
with these wind parameters is Lx ~ lO'^^-lO^^ erg s^^ for binary separations ranging from 
0.5 AU to several AU, assuming a radiative efficiency of 1. 

Currently, there are very few physical examples of SGXBs with X-ray luminosities of 
Lx ~ 10^^ erg s~^ or lower, presumably because X-ray surveys of the Galactic plane are 
biased toward bright SGXBs (e. g., Roche-lobe overflow systems with Lx ~ 10^^ erg s~^; Liu 
et al. 2000). However, a new class of SGXB has recently been unveiled in obscured regions 
of the Galactic plane by hard X-ray surveys, which may provide some insight into X-ray 
production in low-accretion-rate systems: the superfast X-ray transients (SFXT; Negueruela 
et al. 2006a; Sguera et al. 2006) are SGXBs containing hot supergiant primary stars. These 
systems undergo very brief X-ray outbursts, in which the X-ray flux increases by a factor 
of 10^-10^ for durations on the order of hours, in-between much longer states of quiescence 
(Lx ~ 10^^ erg s~^) that last on the order of weeks to a month. Approximately 10 such 
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objects are currently known, many of which contain X-ray pulsars. The outbursts are thought 
to occur when a neutron star in an eccentric orbit nears periastron, and plunges into dense 
inhomogeneities within the hot supergiant's wind. The long quiescent periods ensue as the 
neutron star travels further from the primary; it is this particular state of SFXTs which we 
intend to compare with our sample of Galactic center X-ray sources. For example, during 
their quiescence states the SFXTs XTE J1739-302 and IGR J11215-5952 exhibit hard 
power-law spectra having F 1-2 and X-ray luminosities of Lx ~ (3-4) x 10^^ erg 
within the 0.1-10 keV energy range (Sidoh et al. 2008; Romano et al. 2007, 2009). If these 
particular values of F and Lx for XTE J1739-302 and IGR J11215-5952 were plotted in 
Figure 11, they would be practically indistinguishable from the stars of our sample. Thus, 
the X-ray properties of the Galactic center sources are consistent with the quiescent states 
of SFXTs. Furthermore, quiescent X-ray emission from an SFXT is also consistent with 
Equation (3): IGR J11215— 5952, which contains the Blla star HD 306414, exhibits evidence 
for an orbital period of P ~ 165 days (Romano ct al. 2009). Blla stars have M ^ 40-50 Mq, 
My, ~ jt-\ and Vyj - 1000 km s^^ (Crowther et al. 2006a, 2006b). So, Kepler's law 

implies that this system should have a separation of a ~ 2 AU. Thus, Equation (3) implies 
an X-ray luminosity of Lx ~ 6 x 10^^ erg s~^ for this system, which is a close approximation 
to this source's lowest observed quiescent luminosity of Lx = 4 x 10^^ erg s~^ (Romano et 
al. 2007). 

However, although there has been moderate X-ray variability observed for the several 
sources listed in Table 7, there is no evidence that any of the stars of our sample underwent 
eruptions or outbursts that are on the level known to occur in SFXTs. Thus, we are not 

suggesting that the stars in our sample are SFXTs; rather, we arc examining the X-ray 
properties of quiescent SFXTs to provide some insight into the nature of faint X-ray emission 
from SGXBs with relatively wide binary separations. Indeed, SGXBs with wider and more 
circular orbits could exhibit the quiescent X-ray properties of SFXTs without the periodic 
outbursts. 

In conclusion, the current data are consistent with both the colliding-wind and wind- 
accretion hypotheses. The most effective method of discriminating between these two pos- 
sibilities is to derive a constraint on the system mass functions. Colliding-wind binaries will 
exhibit variable doppler shifts in their spectra on the order of ^ 10-500 km s~^, modulo un- 
favorable binary inclinations or large binary separations, while the doppler signature of the 
gravitational tug of a wind- accreting neutron star (M 2Mq) in the spectrum of a super- 
giant donor star would be more challenging to detect. We plan to perform such experiments 
in the near future. 
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6.2. Spatial Distribution, Environments, and Implications for Steir Formation 

Mid-infrared images from Spitzer /IRAC of the fields containing the massive stellar X- 
ray sources are presented in Figures 13-15. For Figures 13 and 14, three-color images were 
assembled using data from IRAC channels 1 (A3. 6 fj.ni), 3 (A5.8 /xm), and 4 (A8.0 /xm), and 
are displayed with a histogram-equalization contrast stretch to enhance subtle structure in 
surrounding nebulosity. Each image is 144" on a side, which corresponds to a projected 
physical distance of ^5.6 pc, assuming a distance of 8 kpc to the Galactic center. The 
images bear the name of the X-ray source that lies at the field center as well as alphabetical 
labels that show the location of the field in the wide-field A8.0 /im image presented in Figure 
15. 

Several fields in Figures 13-15 contain multiple massive stellar X-ray sources or addi- 
tional, known massive stars: field D of the WC9d star X174519.1 also contains the Ofpe/WN9 
star X174516.1 on the western side of the image and the B2 la star from MMM07 on the east- 
ern side; field F of the 04-61 star XI 74532. 7 (the central source of the HI nebula) contains 
the Of star X174528.6 on the western side of the image and the WN7-8h star X174536.1 
(Mikles et al. 2006) on the eastern side; field G of the O9I-B0I star X174537.3 contains 
the OBe star H5 (Cotera et al. 1999) to the northeast of the field center and the WC8-9 
star WR lOlp (Homeier et al. 2003) to the far northeast; finally, field K of the WC9d star 
X174517.7 also contains the central portion of the Quintuplet cluster toward the north-west 
of the image, which is dominated by the infrared emission of very bright DWCL Quintuplet 
proper members, including X174514.6 (qF231) and X174515.8 (qF211). The LBV Pistol 
Star and its spherical ejection nebula are also apparent in field K. 

Several massive stellar X-ray sources appear to have diffuse mid-infrared nebulae asso- 
ciated with them. In particular, each of the three O supergiants that exhibit the largest A8.0 
fim excess in Figure 8 (X174537.3, X174628.2, and X174528.6) appears to be surrounded by 
such nebulae, which is extended in the case of 174628.2, and more compact in X174537.3, 
and X174528.6. The excess point-source emission may be natal gas and dust in the process 
of being heated and cleared away by the strong radiation fields and winds of the super- 
giants. Although we did not present the SEDs of these stars in Section 5, we did examine 
them, and determined that they appear very similar to the bow-shock source IRS8 near 
the central parsec, which although unresolved by Spitzer, was clearly resolved by Geballe 
et al. (2006) via adaptive-optics imaging. Hence, we speculate that the bright A8.0 //m 
excess from X174537.3, X174628.2 and X174528.6 may be unresolved dust emission from 
bow-shocks or heated dust structures that are similar to that surrounding IRS8. The WN7- 
8h stars X174516.7 and X174712.2 also appear to be associated with bright mid-infrared 
nebulae, although in these cases the nebulae are well resolved from the central stellar point 



-28- 



source. The O9I-B0I star XI 74502. 8 and the WN7-8h star XI 74532. 7 also appear to be 
having strong impact on their surroundings, as evidenced by the bow-shock-hke or shell-hke 
morphology of nearby nebulae. 

The origin of the massive stellar X-ray sources near the Galactic center is currently 
unclear. With the exception of several rather obvious associations with the Arches and 
Quintuplet clusters, it is possible that the relatively isolated stars near the clusters (e. g., 
X174555.3 and X174517.0) originated within them but were dynamically ejected via gravita- 
tional interactions with other massive stars and binaries within the cluster, or were ejected as 
the result of a companion supernova. Indeed, a noteworthy fraction of O stars (~10%-30%) 
in the Galactic disk are known to be runaways, ~10% or more of which are binaries (Gies 
& Bolton 1986; Gies 1987). Runaway Wolf-Rayet-|-0 binaries are also known to exist, such 
as WR 21 and WR 22 in Carina (Moffat et al. 1998). Thus, it is hkely that some massive 
stars and binaries have been dynamically ejected from the Arches and Quintuplet clusters 
over their respective lifetimes of 1-2 Myr (Figer et al. 2002) and 4-6 Myr (Figer, McLean 
&; Morris 1999). It is also tempting to speculate on the ejection scenario for the DWCL 
star X174645.2 (Figure 15, subfield M), owing to its apparent solitude and relatively high 
latitude above the Galactic plane, compared with the other massive stellar X-ray sources. 
Unfortunately, there is no obvious potential point of origin for this evolved star, so its birth- 
place remains elusive. Proper motion measurements performed with adaptive optics, and 
high-resolution spectra, would be required to obtain the kinematic information necessary to 
place meaningful constraints on the origins of the isolated massive stars. 

Alternatively, some of the stars from our sample may be the products of a mode of 
isolated massive star formation, operating in tandem with the formation of dense stellar 
clusters. The relatively loose aggregations of early-0 and WNh stars near the H1-H8 H ll 
regions (Figure 13, fields E and F) and Sagittarius B, and the more evolved WN and WC stars 
southwest of Sagittarius A West might be the products of such a formation mode. However, 
it is also possible that some of these stars are members of stellar clusters less extreme in mass 
and density than the Arches and Quintuplet, which have gone unnoticed owing to confusion 
with the dense stellar background. Indeed, clusters near the Galactic center are subject to 
strong tidal forces, causing them to dissolve on time scales of ~ 10 Myr (Kim et al. 1999), 
while they may fade beyond detectability as a surface density enhancement within only a 
few Myr (Portegies-Zwart et al. 2001). If such hidden clusters exist, narrow-band imaging 
in the vicinity of the newly discovered stars of our sample could betray their presence with 
the detection of additional, coeval emission-line stars. 

Although the origin of the massive stellar X-ray sources reported in this work is puzzling, 
it is clear that significant massive star formation has been occurring in the Central Molecular 
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Zone, outside of the well known clusters, for several Myr, including the formation products 
of the Sagittarius B molecular cloud complex and the newly discovered regions southwest 
of Sagittarius A West. The Chandra observations of the Arches and Quintuplet clusters 
indicate that only a small fraction of massive stars in the Galactic center are detectable as 
hard X-ray sources. Indeed, the Arches and Quintuplet contain over 100 massive stars each 
(M > 2OM0), yet only several (~ 5%) are detected by Chandra as X-ray sources. If roughly 
the same fraction of massive stars is detectable in X-rays throughout the entire Central 
Molecular Zone, we speculate that the existence of ?»20 massive stellar X-ray sources outside 
of the stellar clusters implies the presence of Ri400 additional massive stars that have yet 
to be identified. It is noteworthy that the total inferred Lyman-a photon production rate 
that is required to explain the far-IR luminosity emerging from the central half-kiloparsec 
(~ 10^^ photons s^^) is about twice the amount generated by the three clusters collectively 
(Cox & Laureijs 1989, Figer et al. 2004); the several hundred massive stars implied to exist 
by the results of this work would account for the current deficit. However, dense clusters 
may dynamically produce more hard X-ray sources than less dense associations, in which 
case the number of massive X-ray-emitting stars detected outside of the clusters would imply 
an even larger population of undetected massive stars. 



7. Summciry and Concluding Renicirks 

We have confirmed 16 new massive stars in the Galactic center through infrared spec- 
troscopy of near-infrared counterparts to Chandra X-ray sources. The spectral types include 
early- and late-0 supergiants, as well as WN and WC Wolf-Rayet stars. All of the stars 
exhibit infrared excess, which is attributable to free-free emission from ionized stellar winds, 
supplemented by thermal emission from hot dust in the case of several WC stars. For sev- 
eral cases we modeled the free-free emission and derived stellar wind parameters, producing 
results that are in excellent agreement with previous studies of similar stars. 

The X-ray photometry of the sources exhibit significant fiux at energies above E > 
2 keV, which we regard as an indication of binarity since this property is not typical of 
single massive stars. Our favor currently leans toward the colliding-wind binary hypothesis 
for the majority of our sample, mostly owing to similarities in their X-ray properties with 
known objects from the literature, and the commonality of such objects among massive-star 
associations in the Galaxy. However, the current data arc also consistent with the presence 
of accretion-powered SGXBs in our sample. It should be noted, however, that without 
direct evidence for companions, intrinsic hard X-ray generation from single stars, either via 
magnetically confied winds or some other mechanism, cannot be completely ruled out. Only 



-30- 



a constraint on system mass functions via photometric and radial velocity monitoring will 
allow further discrimination between these alternate hypotheses, and we plan to conduct 
such experiments in the near future. 

The newly discovered Wolf-Rayet stars and O supergiants bring the total number of 
massive stellar X-ray sources to 31 (30 if we regard the X-ray detection of X174502.8 as 
spurious; see Section 3.1), including those within the Arches and Quintuplet clusters. The 

detection of ~ 20 such sources outside of the clusters suggests there may be a large population 
of several hundred massive stars that are currently unidentified in the Galactic center region, 
since wc expect no more than ^5% of massive stars in the Galactic center to be detectable in 
X-rays within our current sensitivity limit. If such a large population exists, it would account 
for the missing Lyman-a photon flux required to explain the far-IR luminosity emerging from 
the central half-kiloparsec. 

Finally, in Paper I we estimated that there are likely to be ^100-300 absorbed X-ray 
sources having real infrared counterparts with Kg < 15.6 mag. Expanding the spectroscopic 
search to include fainter infrared matches to X-ray sources will undoubtedly reveal more 
massive stellar counterparts in the Galactic center, including additional Wolf-Rayet and 
O stars, and HMXBs. Such discoveries will elucidate the origin of these stars and their 
overall distribution throughout the Galaxy's Gentral Molecular Zone. The next generation 
of near-infrared multi-object spectrographs on 8-10 m telescopes, such as FLAMING0S2 on 
Gemini-South (Eikenberry et al. 2008) and MOSFIRE on Keck (McLean et al. 2008), will 
be well suited for this task. 

This research was based upon observations with the Chandra X-ray Observatory. Sup- 
port was provided by the National Science Foundation (grant AST-0406816). J. M. thanks 
Fabrice Martins for providing model stellar spectra of the Arches stars for comparison, and 
comments on the spectral classification of the massive stars presented in this work. We also 
thank Shogo Nishiyama and Tetsuya Nagata for providing near-infrared photometry from 
the SIRIUS survey. 
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Table 1. Infrared Spectroscopic Observations of Counterparts to X-Ray Sources 



Associated X-Ray Source 


Observation Date (UT) 


Telescope/Instrument 


\/S\ 


Seeing 


Sky Cond. 


CXOGC 


J174550, 


,6- 


-285919 


2002 May 8 11:22 


IRTF/SpcX 


2000 




clear 


CXOGC 


J174617, 


,7- 


-285007 


2006 Jul 9 11:06 


AAT/IRIS2 


2400 


1.2" 


clear 


CXOGC 


J174555. 


,3- 


-285126 


2006 Jul 20 08:22 


KeckII/NIRC2 


2200 


0'.'7-0'.'9 


clear 


CXOGC 


J174617. 


,0- 


-285131 


2006 Jul 20 08:48 


KeckII/NIRC2 


2200 


0'.'7-0'.'9 


cleaj: 


CXOGC 


J174656. 


.3- 


-283232 


2007 May 8 13:53 


IRTF/SpcX 


250 


O'.'8-l" 


thin clouds 


CXOGC 


J174711. 


,4- 


-283006 


2007 May 9 11:28 


IRTF/SpcX 


250 


0'.'6-0'.'9 


thin clouds 


CXOGC 


J174703, 


,1- 


-283119 


2007 Jun 15 08:26 


UKIRT/UIST 


2133 


0^4-0^' 5 


thin clouds 


CXOGC 


J174725, 


,3- 


-282523 


2007 Jun 15 09:02 


UKIRT/UIST 


2133 


0"3-0"5 


thin clouds 


CXOGC 


J174628, 


,2- 


-283920 


2007 Jun 16 09:45 


UKIRT/UIST 


2133 


0{'5-0{'6 


clear 


CXOGC 


J174713. 


,0- 


-282709 


2007 Jun 28 06:55 


UKIRT/UIST 


2133 


0'.'9 


thin clouds 


CXOGC 


J174519. 


,1- 


-290321 


2007 Jun 28 07:50 


UKIRT/UIST 


2133 


0'.'9 


thin clouds 


CXOGC 


J174537, 


,3- 


-285354 


2007 Jul 20 06:52 


UKIRT/UIST 


2133 


0('7-0('l 


clear 


CXOGC 


J174532, 


,7- 


-285617 


2007 Jul 22 06:38 


UKIRT/UIST 


2133 


0"7-0"8 


thin clouds 


CXOGC 


J174712, 


,2- 


-283121 


2007 Jul 22 07:36 


UKIRT/UIST 


2133 


0('7-0('8 


thin clouds 


CXOGC 


J174616. 


,6- 


-284909 


2008 May 14 12:51 


AAT/IRIS2 


2400 


0'.'8 


clear 


CXOGC 


J174516. 


.1- 


-290315 


2008 May 16 16:07 


AAT/IRIS2 


2400 


1.3" 


clear 


CXOGC 


J174508. 


,9- 


-291218 


2008 Jun 16 04:01 


SOAR/OSIRIS 


3000 


0'.'8 


30-50% cover 


CXOGC 


J174502. 


,8- 


-290859 


2008 Jun 16 07:29 


SOAR/OSIRIS 


3000 


^.'8 


30-50% cover 


CXOGC 


J174516, 


,7- 


-285824 


2008 Jun 17 05:33 


SOAR/OSIRIS 


3000 


1" 


30-50% cover 


CXOGC 


J174522, 


,6- 


-285844 


2009 Jun 14 01:34 


SOAR/OSIRIS 


1200 


1" 


thin clouds 


CXOGC 


J174645. 


,2- 


-281547 


2009 Aug 05 08:45 


IRTF/SpeX 


1200 


c.'e 


clear 



Table 2. Infrared Photometry of Massive Stellar X-ray Sources 



X-Ray Source 


R.A.IR 


Decl.iR 




J 




H 




Ms. 6 


-M4.B 


Mb. 8 




Ms.o 




(CXOGC J) 




(deg, J2000) 




(mag) 




(mag) 


(mag) 


(mag) 


(mag) 


(mag) 




(mag) 




174502.8- 


-290859'^ 


266, 


,26195 


-29.14986 


13 


.93 ± 0, 


,05 


11.43 ±0.02 


9.88 ±0.01 


8.65 ±0.01 


8.30 ± 0.01 


8.31 ±0.02 








174508.9- 


-291218 


266, 


,28733 


-29.20495 


15 


.12 ± 0, 


,03 


12.61 ±0.02 


11.08 ±0.04 


8.88 ±0.01 


8.42 ± 0.01 


8.15 ±0.02 


8, 


,06 ±0, 


,05 


174516.1- 


-290315 


266, 


,31744 


-29.05430 


11, 


.49 ± 0, 


,02 


9.17 ± 0.03 


7.89 ±0.03 


6.76 ±0.00 


6.28 ±0.01 


6.10 ±0.01 


5, 


,76 ±0, 


,02 


174516.7- 


-285824 


266, 


.31969 


-28.97364 


16, 


.67 ± 0, 


,05 


13.08 ±0.03 


11.09 ±0.02 


9.40 ±0.01 


8.83 ±0.01 


8.68 ± 0.03 


7, 


.38 ± 0, 


,03 


174519.1- 


-290321 


266. 


,32974 


—29.05609 




17.10 




13.30 ± 0.05 


10.40 ± 0.05 


7.64 ± 0.01 


6.78 ± 0.01 


6.55 ±0.01 


6, 


,89 ± 0. 


,03 


174522.6- 


-285844 


266. 


,34453 


—28.97895 








15.21 ±0.06 


12.22 ±0.03 


10.01 ± 0.06 


9.13 ± 0.08 


8.92 ± 0.04 








174528.6- 


-285605 


266, 


.36922 


—28.93476 


14 


.46 ± 0. 


,02 


11.46 ±0.02 


9.72 ± 0.03 


6.58 ±0.01 


5.21 ±0.01 


3.69 ±0.01 


1, 


.43 ± 0, 


,01 


174532.7- 


-285617 


266, 


,38652 


-28.93797 


14, 


.67 ±0, 


,02 


12.13 ±0.01 


10.72 ±0.01 














174536.1- 


-285638 


266, 


,40056 


-28.94405 


15, 


.55 ± 0, 


,03 


12.26 ±0.01 


10.42 ±0.02 


8.82 ±0.01 


8.26 ±0.01 


7.97 ± 0.02 


7, 


.82 ±0, 


,06 


174537.3- 


-285354'' 


266, 


,40538 


-28.89827 


15, 


.75 ± 0. 


,02 


12.80 ±0.01 


11.23 ±0.01 


9.03 ±0.01 


8.18 ±0.01 


6.56 ±0.01 


4, 


,50 ±0. 


,01 


174549.7- 


-284925^^ 


266, 


.45734 


—28.82389 








14.26 ± 0.06 


12.16 ±0.04 














174550.2- 


-284911 


266. 


,45941 


—28.81994 


15, 


.24 ± 0. 


,02 


12.16 ±0.02 


10.47 ±0.01 


9.21 ±0.01 


8.81 ±0.01 


8.82 ± 0.04 








174550.4- 


-284922 


266, 


,46007 


—28.82287 


14, 


.94 ± 0. 


.03 


11.64 ±0.02 


9.87 ± 0.03 














174550.4- 


-284919 


266, 


.46022 


-28.82212 


14 


.48 ± 0. 


,05 


11.45 ±0.02 


9.56 ± 0.06 














174550.6- 


285919 


266, 


,46087 


-28.98878 


15, 


.44 ±0, 


,02 


12.48 ±0.01 


10.85 ±0.01 


9.63 ±0.01 


9.00 ±0.01 


8.63 ±0.02 


7, 


,90 ±0, 


,04 


174555.3- 


-285126 


266, 


,48067 


—28.85738 


15, 


.43 ± 0, 


19 


12.53 ±0.04 


10.97 ±0.07 


9.43 ±0.01 


8.70 ±0.01 


8.54 ±0.03 


8, 


,35 ±0, 


,09 


174614.5- 


-284937 


266, 


,56044 


—28.82688 


14.9 ± 0, 


,2 


12.6 ±0.2 


10.7 ±0.2 














174614.6- 


-284940 


266, 


,56131 


-28.82805 




13.55 




10.69 ±0.11 


7.29 ± 0.02 


3.16 ±0.01 


1.84 ±0.01 


0.81 ±0.01 


0, 


,32 ± 0, 


,01 


174615.1- 


-284932 


266, 


.56318 


-28.82583 


14 


.10 ±0. 


,02 


11.10 ±0.02 


10.59 ±0.20 














174615.8- 


-284945 


266 


,.56611 


-28.829,34 


14, 


.90 ± 0, 


,08 


10.41 ±0.04 


7.24 ±0.03 


3.53 ±0.01 


3.20 ± 0.01 


1.38 ±0.01 


0, 


,98 ± 0. 


,01 


174616.6- 


-284909 


266, 


,56946 


-29.81931 


16 


.53 ± 0, 


,03 


12.99 ±0.01 


11.08 ±0.02 


9.42 ±0.07 












174617.0- 


-285131 


266, 


,57119 


-28.85871 


14, 


.98 ± 0, 


,03 


12.12 ±0.02 


10.49 ± 0.02 


9.01 ±0.01 


8.66 ±0.01 


8.40 ±0.02 


8, 


,50 ±0, 


,05 


174617.7- 


-285007 


266, 


.57430 


-28.83541 


13, 


.99 ±0. 


,05 


10.37 ±0.06 


7.84 ± 0.03 


5.52 ±0.01 


5.02 ±0.01 


3.92 ±0.01 


3, 


.65 ±0, 


,01 


174645.2- 


■281547 


266, 


,68853 


-28.26324 




15.39 




9.96 ± 0.03 


7.18 


saturated 


saturated 


saturated 








174628.2- 


■283920 


266, 


,61793 


-28.65579 


16, 


.99 ± 0. 


,07 


13.36 ± 0.03 


11.49 ± 0.04 














174656.3- 


-283232 


266, 


.73468 


-28.54234 








13.74 ± 0.05 


11.24 ±0.02 


9.50 ±0.01 


8.42 ±0.01 


8.11 ±0.02 


8, 


.25 ±0, 


,05 


174703.1- 


-283119 


266, 


,76307 


-28.52225 


16, 


.23 ± 0, 


,03 


13.03 ±0.01 


11.27 ±0.01 


7.52 ±0.01 


6.31 ±0.01 


5.00 ±0.01 


3, 


,45 ± 0, 


,01 


174711.4- 


283006 


266, 


,79782 


-28.50194 


16, 


.56 ± 0, 


,06 


12.72 ±0.02 


10.54 ±0.01 


9.29 ±0.01 


8.53 ±0.01 


8.33 ±0.02 


8, 


.51 ±0, 


,06 


174712.2- 


-283121 


266, 


,80104 


-28.52267 


17, 


.06 ± 0. 


,07 


13.07 ±0.01 


10.78 ±0.02 


9.01 ±0.02 


8.21 ±0.01 


7.42 ±0.03 








174713.0- 


-282709 


266, 


.80423 


-28.45251 








14.22 ±0.01 


11.86 ±0.01 


10.57 ±0.01 


9.65 ±0.01 


9.34 ± 0.03 








174725.3- 


-282523 


266, 


,85567 


-28.42304 








13.37 ±0.01 


11.30 ± 0.04 


10.24 ±0.01 


9.59 ± 0.01 











Note. — All known massive stars in the Galactic center with X-ray counterparts are listed, in addition to those discovered in this work. IR 
source positions and JHKg photometry are taken from Paper I. The A3.6-A8.0 /xm photometry is from the Spitzer /IKAC point-source catalog 



of Ramirez et al. (2008). Uncertainties are included where available. 

*X-ray source considered spurious and not included in main X-ray catalog of Muno et al. (2009). This may be an X-ray source near our 
limit of sensitivity (see text). 

''This source did not make the final candidate counterpart list in Paper I; it has an IR positional error greater than 0.1". We suspect this 
large error is a result of hot dust emission contribution to the Ks-hand, since this source also has a very bright A8 /im counterpart. 

'^JHK photometry for this source from Figer et al. (1999). Note that in this case a /sT-band measurement is presented, rather than the Ks 
band. 

''as /im data from the MSX Point Source Catalog, assuming the zero-point flux of 58.4 Jy from Cohen et al. (2000). 
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Table 3. Spectral Types and Absolute Photometry for WN and O Supergiant X-ray Sources 



X-Ray Source 


Spectral 


Ks 


{J-Ks)o 




J -Ks 












BCjf 


Log 


(CXOGC J) 


TvDe 


V **&/ 












rmaff^ 


( maff^ 


^^maEF^ 


^^maEF^ 




174502.8-290859 


O9I-B0I 


9.88 


-0.21 


-0.10 


4.05 


1.55 


2.10 


2.38 


2.24 


-6.86 


-3.7 


6.12 


174508.9-291218 


WN7 


11.08 


0.36 


0.26 


4.04 


1.53 


1.82 


1.83 


1.83 


-5.25 


-3.5 


5.40 


174516.1-284909 


Ofpe/WN9 


7.89 


0.13 


0.11 


3.60 


1.28 


1.71 


1.68 


1.70 


-8.31 


-2.9 


6.38 


174516.7-285824 


WN7-8h 


11.09 


0.13 


0.11 


5.58 


1.99 


2.69 


2.71 


2.70 


-6.11 


-4.3 


6.06 


174522.6-285844 


WN5-6b 


12.22 


0.36 


0.26 




2.99 




3.93 


3.93 


-6.21 


-3.8 


5.90 


174528.6-285605 


WN8-9h 


9.72 


0.13 


0.11 


4.74 


1.74 


2.28 


2.35 


2.31 


-7.09 


-4.3 


6.46 


174532.7-285126 


04-6I 


10.72 


-0.21 


-0.10 


3.95 


1.41 


2.06 


2.17 


2.11 


-5.89 


-4.4 


6.02 


174536.1-285638 


WN8-9h 


10.42 


0.13 


0.11 


5.13 


1.84 


2.47 


2.49 


2.48 


-6.56 


-4.3 


6.24 


174537.3-285354 


O9-B0I 


11.23 


-0.21 


-0.10 


4.52 


1.57 


2.34 


2.40 


2.37 


-5.64 


-3.7 


5.64 


174549.7-284925 


OB 


12.16 


-0.21 


-0.10 




2.10 




3.17 


3.17 


-5.51 


-3.7 


5.58 


174550.2-284911 


WN8-9h 


10.47 


0.13 


0.11 


4.77 


1.69 


2.29 


2.28 


2.29 


-6.32 


-4.3 


6.15 


174550.4-284922 


WN8-9h 


9.87 


0.13 


0.11 


5.07 


1.77 


2.44 


2.39 


2.42 


-7.05 


-4.3 


6.44 


174550.4-284919 


WN8-9h 


9.56 


0.13 


0.11 


4.92 


1.89 


2.37 


2.56 


2.46 


-7.40 


-4.3 


6.58 


174550.6-285617 


WN7 


10.85 


0.36 


0.26 


4.59 


1.63 


2.09 


1.97 


2.03 


-5.68 


-3.5 


5.57 


174555.3-285126 


WN5-6b 


10.97 


0.36 


0.26 


4.46 


1.56 


2.03 


1.87 


1.95 


-5.48 


-3.8 


5.61 


174614.5-284937 


OBI 


10.7 


-0.21 


-0.10 


4.2 


1.9 


2.18 


2.88 


2.53 


-6.33 


-3.7 


5.91 


174615.1-284932 


OBI 


10.59 


-0.21 


-0.10 


3.51 


0.51 


1.84 


0.88 


1.36 


-5.27 


-4.3 


5.73 


174616.6-284909 


08-9I 


11.08 


-0.21 


-0.10 


5.45 


1.91 


2.80 


2.89 


2.85 


-6.27 


-3.7 


5.89 


174617.0-285131 


06If+ 


10.49 


-0.21 


-0.10 


4.49 


1.63 


2.32 


2.49 


2.41 


-6.42 


-4.3 


6.19 


174628.2-28.3920 


04-61 


11.49 


-0.21 


-0.10 


5.50 


1.87 


2.82 


2.84 


2.83 


-5.84 


-4.4 


6.00 


174656.3-283232 


WN8-9h 


11.24 


0.13 


0.11 


0.00 


2.50 


0.00 


3.44 


3.44 


-6.70 


-4.3 


6.30 


174703.1-283119 


04-6I 


11.27 


-0.21 


-0.10 


4.96 


1.76 


2.55 


2.68 


2.62 


-5.85 


-4.4 


6.00 


174711.4-283006 


WN8-9h 


10.54 


0.13 


0.11 


6.02 


2.18 


2.91 


2.98 


2.95 


-6.91 


-4.3 


6.38 


174712.2-283121 


WN7-8h 


10.78 


0.13 


0.11 


6.28 


2.29 


3.04 


3.14 


3.09 


-6.81 


-4.3 


6.34 


174713.0-282709 


WN7-8h 


11.86 


0.13 


0.11 


0.00 


2.36 


0.00 


3.24 


3.24 


-5.88 


-4.3 


5.97 


174725.3-282709 


04-6I 


11.30 


-0.21 


-0.10 


0.00 


2.07 


0.00 


3.12 


3.12 


-6.32 


-4.4 


6.19 



Note. — Intrinsic colors and BCk values were adopted from Crowther et al. (2006a) for all WN X-ray sources. For the 04— 61 stars, intrinsic 

colors and BC/f values were adopted from Martins & Plez (2006) by averaging the values for spectral type 04— 61, and adopting the values for 
an 09.5I star to represent that of the 09I— BOI stars. We did not include WCd or DWCL stars in this analysis, since we do not have knowledge 
of their intrinsic stellar colors. 
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Table 4. Free-free Model Parameters 



Star 


Spectral 




R* 


(3.6, 4.5, 5.8, 8 


1.0 fim) 




ATe X 10" 


% 


At=1 




Type 


(kK) 


(Re) 


(Jy) 




(K) 


(cm~^) 


(Rq) 




X174516.1 


Ofpe/WN9 


20 


79 


1.60 : 1.74 : 1.87 


: 2.00 


8100 


0.8 


246 


5.2 


X174656.3 


WN8-9h 


30 


33 


0.29 : 0.26 : 0.21 


: 0.14 


7000 


1.3 


131 


4.7 


X174712.2 


WN7-8h 


35 


24 


0.31 : 0.31 : 0.26 


; 0.18 


8500 


1.1 


145 


5.9 


X174555.3 


WN5-6b 


40 


12 


0.16 : 0.16 : 0.14 


: 0.10 


11000 


1.5 


94 


6.4 



Table 5. Chandra X-ray Data for Massive Stellar X-ray Sources in the Galactic Center 



X-Ray Source 


R. 


Ax 


Decl.x 


''X 


Exposure 


Cnet'* 


/~< hard 
•^net 


not (0.5 


-8 


keV) 




HRO 






HR2 




(E) 


Logix 


Spectral 


Ref.» 


Ref.'' 


(CXOGC J) 




(deg, J2000) 


(arcsec) 


(sec) 






(cm" 


-2 


s~ 
















(keV) 


(erg s 




Type 


X 


IR 


(1) 




(2) 


(3) 


(4) 


(5) 


(6) 


(7) 






(8) 






(9) 






(10) 




(11) 


(12) 


(13) 


(14) 


(15) 


174608.9-291218 


266, 


.28731 


-29, 


.20474 


0.87 


96809 


< 5.8 




3, 


.39 


X 


10 


-6 


1 


.oo:« 


.00 
.95 


0, 


.37+4 


.24 
.24 


6.3 


32, 


.9 


WN7 


1 


1 


174616.1-290315 


266, 


.31750 


-29 


.05430 


0.44 


936670 




748.oti:i 


3, 


.42 


X 


10 


-6 


0, 




07 
.07 


-0, 


.33tS 


. 10 
.10 


3.3 


32, 


.0 


Ofpe/WN9 


2 


2 


174616.7-285824 


266, 


.31959 


-28, 


.97368 


0.66 


1011840 


< 7.1 


660l22,o 


3, 


.56 


X 


10 


-7 


1 




.00 

,59 


-0, 


.26+4 


,40 
.44 


4.0 


32, 


,0 


WN7-8h 


1 


1 


174519.1-290321 
174522.6-285844 


266 
266, 


,32983 
.34457 


-29 
-28, 


.05605 
.97898 


0.51 
0.63 


1011810 
1011840 


< 14.3 

< 6.3 


197 5+26.5 

4i.3t;i- 


1 
2 


.05 
.30 


X 
X 


10 
10 


-6 
-7 


1 
1 


.oo:» 


00 
,37 
,00 
,95 




-0, 


.16+4 
.33tS 


.17 
.16 
.34 
.44 


5.0 
4.4 


32 
31, 


4 
,8 


WC9d 
WN5-6b 


1 
1 


1 
1 


174628.6-286606 


266, 


.36926 


-28, 


.93483 


0.48 


1011840 


8.3+_l,l 


178.3l|;«« 


9, 


.00 


X 


10 


-7 


0, 


.82l°- 


12 
,13 


-0, 


■32+4 


,18 
.18 


3.7 


32, 


.4 


Of 


2 


6, 8 


174632.7-286617 


266 


.38650 


-28, 


.93807 


0.60 


1011840 




158.5+25;g 


7, 


.18 


X 


10 


-7 


0, 




22 
.16 


-0, 


.13+4 


,21 
.19 


3.9 


32. 


.3 


04-6I 


1 


1 


174536.1-285638 


266, 


.40059 


-28, 


.94407 


0.37 


1011830 


216.7^2=;° 


4860. Ot^l" 


1 


.97 


X 


10 


-5 


0, 


.77+4 


03 

,02 


-0, 


.20+4 


,02 
.02 


3.9 


33, 


.7 


WNS-9h 


3 


3 


174537.3-285354 


266 


.40558 


-28, 


,89844 


0.74 


1011830 




91 0+^*'" 


4, 


.17 


X 


10 


-7 


0, 


.5lt- 


29 
31 


-1, 


.00+4 


,64 
.00 


3.4 


32, 


,0 


O9I-B0I 


1 


1 


174549.7-284925 


266, 


.45722 


-28, 


.82385 


0.41 


110174 


< 3.6 


42 1+130 


1 


.62 


X 


10 


-6 


1 


.001° 


.00 
.66 


0, 


•llt°o 


,31 
.32 


4.9 


32, 


.6 


OB 


4 


6 


174660.2-284911 


266, 


.45939 


-28, 


.81994 


0.33 


110174 


25.7+4:1 


649.lt|^;^, 


3, 


.00 


X 


10 


-5 


0, 


77+0. 


07 

.07 


-0, 


,02±« 


.08 
.08 


4.1 


33 


.9 


WN8-9h 


4 


6, 8, 9 


174660.4-284922 


266, 


.46004 


-28, 


.82285 


0.33 


110174 


43 4+11" 


958.0tt»;S 


3, 


.97 


X 


10 


-5 


0, 


77+° ■ 
■ ' ' -0 


05 

,06 


-0, 


lj + 


.07 
.07 


4.0 


34, 


.0 


WN8-9h 


4 


6, 8, 9 


174650.4-284919 


266, 


.46025 


-28, 


.82210 


0.33 


110174 


c+9.8 


681.4t|l;| 


2, 


.29 


X 


10 


-5 


0, 


7g+0, 


06 

,06 


-0 


27+" 
.2.1 _Q 


.09 
.09 


3.7 


33, 


,8 


WNS-9h 


4 


6, 8, 9 


174550.6-285919 


266 


.46088 


-28 


,98879 


0.43 


1011850 




43 7+16.0 
' -21.7 


2 


.06 


X 


10 


-7 


-1 


.oo;i- 


2 1 
00 


-0 




.46 
.58 


4.0 


31 


,7 


WN7 


2 


8 


174555.3-285126 


266, 


.48075 


-28, 


.85738 


0.49 


208788 


< 5.6 


86.2ll8| 


1 


.65 


X 


10 


-6 


1 


.oo-J> 


,00 
,30 


-0, 


.6^+4 


.25 
.28 


3.5 


32, 


,6 


WN5-6b 


1 


1 


174614.6-284937 


266, 


.56058 


-28, 


.82698 


0.60 


168922 


46 1+120 


20.7±|-« 


1 


.86 


X 


10 


-6 


-0, 


.60+4 


18 
.18 


-1, 


.oo±° 


.90 
.00 


2.1 


32, 


.7 


OBI 


4 


11 


174614.6-284940 


266, 


.56113 


-28, 


.82794 


0.62 


168922 


< 3.4 


27 9+^ '^ 
^'■"-9.7 


8, 


.90 


X 


10 


-7 


1 


.00-4 


,00 
,32 


-1 


.00+4^ 


.82 
.00 


3.3 


32, 


.3 


DWCL 


4 


7 


174615.1-284932 


266, 


.56311 


-28, 


.82578 


0.61 


168922 


< 6.9 


61 4+l'''l 


1 


.36 


X 


10 


-6 


1 


.00-4 


,00 
,23 


-0 


.r,i+4 


.35 
.36 


3.7 


32, 


.6 


OBI 


4 


7 


174615.8-284945 


266 


.56610 


-28 


,82931 


0.49 


158922 


< 4.5 




1 


,80 


X 


10 


~6 


1 


.00-4 


,00 
,73 


-0, 


.22+4 


.23 
.23 


4.5 


32, 


,7 


DWCL 


4 


7 


174616.6-284909 


266, 


.56939 


-28, 


.81920 


0.75 


208347 


< 6.0 


31 1 + "-^ 


4, 


.91 


X 


10 


-7 


1 


.00-4 


,00 
,48 


-9, 


.00-4 


.00 
.00 


2.9 


32, 


.1 


08-9I 


4 


7 


174617.0-286131 


266 


.57102 


-28, 


.85867 


0.66 


208347 


< 8.4 


114.5t22-5 


2, 


.38 


X 


10 


-6 


1 


.00-4 


,00 
.22 


-0, 


27+0 
■'"-0 


.25 
.26 


3.9 


32, 


.8 


06If+ 


1 


1 


174617.7-286007 


266, 


.57417 


-28, 


.83537 


0.60 


168922 


< 4.6 




1 


.60 


X 


10 




1 


■ooio 


.45 


-0 


48 + " 


.27 


3.6 


32, 


.6 


WC9d 


1 


1 


174628.2-283920 


266, 


.61791 


-28, 


.65569 


0.61 


49481 


< 2.5 


24.ol|;i 


3, 


.62 


X 


10 


-6 


1 


.00-4 


,00 
,52 





12 + " 


.41 
.41 


4.4 


33, 


,0 


04-6I 


1 


1 


174645.2- 281547 

174656.3- 283232 


266 
266, 


.68842 
.73465 


-28, 
-28, 


.26319 
.54230 


1.04 
0.67 


31254 
157384 


< 3.5 

< 4.3 


96^ 1 +27-9 
135.8_2Q g 


6 
3, 


,83 
.91 


X 
X 


10 
10 


-5 
-6 


1 
1 


■Ot'-O 

.00-4 


,00 
,38 
,00 
,41 




-0 


.2S+4 
.09+4 


.10 
.10 
.17 
.17 


5.2 
4.4 


34, 
33, 


,2 
,0 


DWCL 
WN8-9h 


1 


5 
1 


174703.1-283119 


266 


.76298 


-28, 


.52222 


0.7 


107903 


< 2.4 


42 4+11''' 


1 


.83 


X 


10 


-6 


1 


.00-4 


,00 
,45 


-0 


16 + " 


.31 
.30 


4.2 


32, 


.7 


04-6I 


1 


1 


174711.4-283006 
174712.2-283121 


266, 
266, 


.79763 
.80097 


-28, 
-28, 


.50194 
.52258 


0.60 
0.71 


118333 
118333 


< 2.37 

< 2.5 


36.12+llg3 
45 3+12 6 


1 
1 


.13 
.74 


X 

X 


10 
10 


-6 
-6 


1 
1 


.00-4 
.00-4 


,00 
,39 
,00 
,66 


-0, 
-0, 


19+" 


.35 
.32 
.29 
.32 


3.6 
4.6 


32, 
32, 


.4 
.6 


WN8-9h 
WN7-8h 


1 
1 


1 
1 


174713.0-282709 


266 


.80425 


-28 


,45251 


0.45 


107907 


< 1.7 


17 7+''" 


6, 


,77 


X 


10 


-7 


1 


.00-9 


,00 
,19 


-0, 


,27i» 


.39 
.48 


4.3 


32, 


,2 


WN7-8h 


1 


1 


174725.3-282523 


266, 


.85559 


-28, 


.42301 


0.60 


128753 


< 3.5 




2 


.21 


X 


10 


-6 


1 




.00 
.92 


-0, 


,02+"„ 


,25 
.24 


4.6 


32, 


.7 


04-6I 


1 


1 



Note. — The tabic columns contain the following: 

(I) The source name, which was derived from the coordinates of the sources based on the lAU format, in which least-significant figures are truncated (as opposed to rounded). The names 
should not be used as the locations of the sources. 

(2—3) The right ascension and declination of the X-ray source, in degrees (J2000). 

(4) The 95%-confidence uncertainty in the position (the error circle). 

(5) The total exposure time in seconds. 

(6) The net number of counts in the soft U.5— 2.U keV energy band, with 90%-confidence uncertainties. 

(7) The net numbers of counts in the hard 2.0—8.0 keV energy band, with 90%-confidence uncertainties. 

(8) The total 0.5—8.0 keV flux in units of photons cm ^ 

(9) The soft color, defined as HRO = (h — s)/{h -\- s), where h and s is the flux in the 2.0-3.3 keV and 0.5-2.0 keV bands, respectively. The upper and lower 90%-confidence uncertainty is 
included. An HR value of —9.00 indicates that the sources was not detected in either of the h or s bands. HR uncertainties of —9.00 are considered null, and are present whenever a source 
is detected in only one band. 

(10) The hard color, defined as HR2 = {h' — s')/ (h' -\- s'), where h' and s' is the flux in the 4.7-8.0 keV and 3.3-4.7 keV bands, respectively. The upper and lower 90%-confldence uncertainty 
is included. 

(II) The average energy per photon in units of keV. This value can be multiplied by the photon flux in column 8 for an estimate of energy flux. 



(12) The estimated X-ray luminosity of the source (0.5—8.0 keV), assuming A^j^ — 6 X 1(J^^ cm~^ and kT — 2 keV. Owing to our insensitivity to soft X-ray photons between 0.5- 
the uncertainties are significant, so, these values should be met with caution (see text). (13) The spectral type of the stellar counterpart. 

(14) References to the first published association of a massive star with an X-ray source. 

(15) References to the spectral classification of a source. 

^The following authors first recognized a massive-star as an X-ray source: (1) This work or Mauerhan et al. 2007; (2) Muno et al. 2006b; (3) Mikles et al. 2006; (4) Yusef-Zad 
2002; (5) Hyodo et al. 2008. 

^The following authors determined the spectral type of the associated massive star from infrared observations: in addition to the authors listed in a), (6) Figer et al. 2002; (7) 
al. 1999; (8) Cotera et al. 1996 or Cotera et al. 1999; (9) Martins et al. 2008; (10) Clark et al. 2009; (11) Liermann et al. 2009. 



I 
I 
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Tabic 6. Lx Compared with Other Authors 



Source 
(CXOGC J) 


kT (us) 
(keV) 


A^H (us) 
(cm~2) 


LogLx (us) 
(erg s-i) 


kT (them) 
(keV) 


Nu (them) 
(cm~2) 


LogLx (them) 
(erg s-l) 


ALx 
(dex) 


Ref. 


174528.6-285605 
174550.2-284911 
174550.4-284922 
174550.4-284919 
174550.6-285919 
174645.2-281547 


2 
2 
2 
2 
2 
2 


6 
6 
6 
6 
6 
6 


32.4 
33.9 
34.0 
33.8 
31.7 
34.2 


1 9+0.2 

o r,c+2.62 

o.zo_-^ 24 

-,+0.58 
^•^-0.34 

1 07+0.39 
-'-■°'-0.32 

2 

0+0.6 


8.lll-i 
7.3tl;? 
6 

2313 


33.1 
33.7 
33.9 
34.0 
31.7 
34.5* 


0.7 
0.2 
0.1 
0.2 
0.0 
0.3 


2 
3 
3 
3 
2 
1 



Note. — This table compares the X-ray luminosities derived via (1) X-ray photometry, adopting constant 
average values for plasma temperature (kT) and interstellar absorption (Afjj), and (2) X-ray spectral fitting with 
temperature and absorption as free parameters. Only non-variable X-ray sources were used for this comparison. 
References: (1) Hyodo ct al. 2008; (2) Muno et al. 2006; (3) Wang et al. 2006. 

*The luminosity of this source is for the 2.0—8.0 keV band. 



Table 7. Massive Stars Exhibiting X-Ray Variability 



Record 
No. 


Source 
(CXOGC J) 


Min. ObsID 


UT Date 


F ■ 

^ min 

(10-6 cm-2 s-l) 


Max. ObsID 


UT Date 


-Pmax 

(10-6 cm-2 s-l) 


^max/^min 


247 
1733 
8598 
7624 


174516.1-290315 

174536.1- 285638 

174712.2- 283121 
174614.5-284937 


6113 
2952 
944 
945 


2005 Feb 27 06;26;24 
2002 Mar 23 12:25:04 
2000 Mar 29 09:44:36 
2000 Jul 07 19:05:19 


<2.8 
7.1 ± 1.9 
1.1 ±0.3 
1.7 ±0.9 


4684 
5950 
2277 
4500 


2004 Jul 06 22:29:57 

2005 Jul 24 19:58:27 

2001 Jul 16 11:52:55 
2004 Jun 09 08:50:32 


9.4 ± 1.0 
24.0 ± 1.5 
4.1 ± 1.9 
2.7 ±0.4 


o (,-1.0 

■^•"1.8 
^■°0.6 



Note. — The columns give: the record number from the catalog of Muno et al. (2009), the source name, the ObsID and UT date in which the minimum 
flux was observed, the minimum flux, the ObsID and UT date in which the maximum flux was observed, and the maximum flux. The final column contains 
the ratios of the maximum and minimum observed fluxes, with upper and lower uncertainties. Uncertainties of —1.0 imply that the flux ratio is a lower limit. 
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Fig. 1. — K-hand spectra of late-O/early-B supergiant X-ray sources in the Galactic center, 
exhibiting absorption hnes of H I and He I, and a weak emission feature from N ill near A2.11 
/xm. The spectra of X174537.3 and X174502.8 resemble that of O9I-B0I stars, while the 
presence of He II absorption and a relatively shallow Br7 line in the spectrum of X174616.6 
are indicative of a slightly earlier spectral type in the range of 08-91. 
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Fig. 2. — X-band spectra of early 04-61 supergiant X-ray sources, dominated by emission 
from the complex of He I, N ill, C ill, and O ill near A2. 112-2. 115 //m. Weak features of 
Br7 emission or absorption and He II absorption at A2.19 are also present. A model 
spectrum of an 04-61 star is plotted near the bottom of the figure for comparison (Martins 
et al. 2008). 
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Fig. 3.— iT-band spectra of WN8-9h X-ray sources X174656.3 and X174711.4, and the 
061f"^ star X174617.0 from Mauerhan, Muno, & Morris (2007). Strong Br7 emission is 
accompanied by the emission complex of He I, N III, C III and O III near A2. 112-2. 115 fim. 
A model spectrum of a WN8-9h star is plotted near the bottom of the figure for comparison 
(Martins et al. 2008). 
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Fig. 4. — i^-band Spectra of WN7-8h X-ray sources, which exhibit stronger He ll emission 
at A2.189 /im, relative to the WN8-9h stars presented in Figure 3. 
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Fig. 5. — i^-band spectra of hydrogen-deficient Wolf-Rayet X-ray sources. Broad emission 
lines of He I and He ll dominate the spectra, indicating high temperatures and fast extended 
winds. X174550.6 and X174508.9 have spectral types of WN7; X174555.3 and X174522.6 
have spectral types in the range of WN5-6b. 
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Fig. 6. — i^-band spectra of WC counterparts to Galactic center X-ray sources. The weak- 
ness of the emission lines is attributable to dilution by strong excess continuum emission 
from hot dust, and perhaps from a supergiant companion. The spectrum of X174645.2 is an 
extreme example of dust dilution, and is completely featureless within the noise level. 
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Fig. 7. — X-band spectra of the counterpart to X-ray source CXOGC J174516. 1—281531. 
Strong emission lines of Br7 and He I P-Cygni at A2.058 /xm, are accompanied by weaker 
lines of Mg ll near A2.14 /xm, and weak He I and He ll P-Cygni profiles near A2.11 /xm and 
A2.19 /xm, respectively. The spectrum resembles that of Ofpe/WN9 stars. 
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Fig. 8. — Infrared color-color diagram of massive stellar X-ray sources near the Galactic 
center, plotted with la error bars. Almost all of the massive stars (large symbols) exhibit 
significant long-wavelength excess compared to field stars from the SIRIUS near-infrared 
survey (small dots). 
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Fig. 9. — Extinction-corrected Al-10 fim photometry (filled squares) and model SEDs for 
selected WN stars exhibiting free-free emission in the infrared. For the model SEDs, the 
dashed lines represent the Rayleigh-Jeans portion of the stellar continuum; the dash-dotted 
lines represent the free-free emission spectrum; and the solid lines are the model sum. 
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Fig. 10. — Extinction-corrected Al-10 /im SEDs for new WC9d stars X174519.1 and 
X174617.7 (filled triangles), compared with the dust-enshrouded, extreme DWCL star qF231 
of the Quintuplet cluster (squares), the single WC9 star HDM13 and the WC9d star WR 
59 (circles). Rayleigh-Jeans curves (dotted lines) appropriate for WC9 stars {T^ff = 40 kK) 
are scaled to the J-band data points for comparison with the SED curves. The SEDs of 
X174519.1 and X174617.7 exhibit an infrared excess similar to that of WR 59, consistent 
with the WC9d class. 
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Fig. 11. — Photon flux (0.5-8.0 keV) vs. hard X-ray color (HR2) for confirmed massive-star 
counterparts to Chandra X-ray sources. Sources not detected at energies above 3.3 keV 
were assigned HR2=— 1.1 for the figure, but have HR2=— 9.00 in Table 5. The massive 
stars (large dots) are systematically brighter and softer in X-rays than the field population 
of X-ray sources (small grey points). Sources associated with the Arches (pentacles within 
circles) and Quintuplet (pentacles) clusters are marked separately for comparison. The data 
are accompanied by models of optically thin thermal plasma (solid lines) and non-thermal 
power-law emission (dotted lines) absorbed by a hydrogen column of A^h = 6.0 x 10^^ cm~^. 
The implied X-ray luminosities of the massive stars are in the range Lx ~ 10^^-10'^^ erg s~^ 
(0.5-8.0 keV). 
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Fig. 12. — Scatter plots illustrating the relationship between the X-ray luminosity (0.5-8.0 
keV) and the stellar bolometric luminosity of the massive stellar X-ray sources near the 
Galactic center. A hnear fit to the data (dashed line) yields a value of Lx/ L^oi ~ lO~^-^='=°-^, 
which is consistent with the canonical value of 10~^ (sohd line). 
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Fig. 13. — S])itzeT j'^KKQ, images of fields containing massive stellar X-ray sources. Red, 
green, and blue correspond to 8.0, 5.8, and 3.6 [im. images, respectively. Each image spans 
144" on a side, and is centered on the labeled sources. Other known massive stars, if present 
within these fields, are marked by perpendicular pointers. In all images, north is up and 
east is to the left. Image contrast was scaled by histogram equalization to help reveal faint 
structure, so the colors of the images are not necessarily illustrative of source SEDs. 
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Fig. 14. — Spitzer /IRKC images of fields containing additional massive stellar X-ray sources. 




Fig. 15. — Spitzer /IRKC A8 /xm image of the Galactic center region (Stolovy et al. 2006; Arendt et al. 2008), oriented 
so the Galactic plane is horizontal. The field of each massive stellar X-ray source studied in this work is marked, 
matching the individual, zoomed color images in Figures 13 and 14. 



